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Abstract We review the stability properties of several discretizations of the Helmholtz 
equation at large wavenumbers. For a model problem in a polygon, a complete k- 
explicit stability (including A:-explicit stability of the continuous problem) and con- 
vergence theory for high order finite element methods is developed. In particular, 
quasi-optimality is shown for a fixed number of degrees of freedom per wavelength 
if the mesh size h and the approximation order p are selected such that kli / p is suf- 
ficiently small and p = 0{logk), and, additionally, appropriate mesh refinement is 
used near the vertices. We also review the stability properties of two classes of nu- 
merical schemes that use piecewise solutions of the homogeneous Helmholtz equa- 
tion, namely. Least Squares methods and Discontinuous Galerkin (DG) methods. 
The latter includes the Ultra Weak Variational Formulation. 



1 Introduction 

A fundamental equation describing acoustic or electromagnetic phenomena is the 
time-dependent wave equation 



w 



d 



given here for homogeneous, isotropic media whose propagation speed of waves 
is c. It arises in many applications, for example, radar/sonar detection, noise fil- 
tering, optical fiber design, medical imaging and seismic analysis. A commonly 
encountered setting is the time-harmonic case, in which the solution w (and corre- 
spondingly the right-hand side g) is assumed to be of the form Re (e^"'"M(x)) for 
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a frequency co. Upon introducing the wavenumber k = (o/c and the wave length 
X :~ 2n/k, the resuhing equation for the function u, which depends solely on the 
spatial variable x, is then the Helmholtz equation 

-Au-k^u = f. (1) 

In many high frequency situations of large k the solution u is highly oscillatory 
but has some multiscale character that can be captured, for example, by means of 
asymptotic analysis; a classical reference in this direction is Q. 

In this article, we concentrate on numerical schemes for the Helmholtz equation 
at large wavenumbers k. Standard discretizations face several challenges, notably: 

(I) For large wavenumber k, the solution u is highly oscillatory. Its resolution, 
therefore, requires fine meshes, namely, at least = degrees of freedom, 
where d is the spatial dimension. 
(II) The standard //'-conforming variational formulation is indefinite, and stabil- 
ity on the discrete level is therefore an issue. A manifestation of this problem 
is the so-called "pollution", which expresses the observation that much more 
stringent conditions on the discretization have to be met than the minimal 
= 0[¥') to achieve a given accuracy. 

The second point, which will be the focus of the article, is best seen in the following, 
one-dimension example: 

Example 1.1. For the boundary value problem 

-u"-k^u^\ in (0,1), h(0)=0, M'(l)-iyt«(l) =0, (2) 

we consider the /i-version finite element method (FEM) on uniform meshes with 
mesh size h for different approximation orders p G {1,2,3,4} and wavenumbers 
k € {1,10,100}. Fig. [U shows the relative error in the //'(i2)-semi norm (i.e., 
|w — MA'|//i(i2)/|M|//i(^2)' where un is the FEM approximation) versus the number 
of degrees of freedom per wavelength A^;^ ■ = N/X = 2nN/k with being the di- 
mension of the finite element space employed. We observe several effects in Fig.[T] 
Firstly, since the solution « of ^ is smooth, higher order methods lead to higher 
accuracy for a given number of degrees of freedom per wavelength than lower order 
methods. Secondly, asymptotically, the FEM is quasioptimal with the finite element 
error \u — MA'|//i(i2) satisfying 

for a constant Cp independent of k. Thirdly, the performance of the FEM as mea- 
sured in "error vs. number of degrees of freedom per wavelength" does depend on 
k: As k increases, the preasymptotic range with reduced FEM performance becomes 
larger Fourthly, higher order methods are less sensitive to k than lower order ones, 
i.e., for given k, high order methods enter the asymptotic regime in which (|3) holds 
for smaller values of Nx than lower order methods. ■ 




The behavior of the FEM in Example I l.ll has been analyzed in 14711491 . where error 
bounds of the form (see ||47l Thm. 4.27]) 

|«-«iv|//.(i2) <C;,(l+F+V)/z'>U.+,(^) (4) 

are established for a constant Cp depending only on the approximation order p. In 
this particular example, it is also easy to see that \u\ni>+\n)l\'AH^{n) ^ ^o that 
(lUl can be recast in the form 

This estimate goes a long way to explain the above observations. The presence of 
the factor 1 + ^A^^ explains the "pollution effect", i.e., the observation that for 
fixed Nx, the (relative) error of the FEM as compared with the best approximation 
(which is essentially proportional to A^^ in this example) increases with k. The es- 
timate (|5]l also indicates that the asymptotic convergence behavior (|3]l is reached for 
Ni = 0{k^lP). This confirms the observation made above that higher order meth- 
ods are less prone to pollution than lower order methods. Although Example 11.11 
is restricted to ID, similar observations have been made in the literature also for 
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multi-d situations as early as [[T5| . We emphasize that for uniform meshes (as in Ex- 
ample [TTT) or, more generally, translation invariant meshes, complete and detailed 
dispersion analyses are available in an /i-version setting, ||2]|27]|42l|49l, and in a 
p/Zip-setting, [2-4 1, that give strong mathematical evidence for the superior ability 
of high order methods to cope with the pollution effect. 

The present paper, which discusses and generalizes the work II61II62I . proves that 
also on unstructured meshes, high order methods are less prone to pollution. More 
precisely, for a large class of Helmholtz problems, stability and quasi-optimality is 
given under the scale resolution condition 

kh 

— <c\ and p>C2\ogk, (6) 
P 

where c\ is sufficiently small and C2 sufficiently large. Forpiecewise smooth geome- 
tries (e.g., polygons), additionally appropriate mesh refinement near the singularities 
is required. 

We close our discussion of Example 1 1 . 1 1 bv remarking that its analysis and, in 
fact, the analysis of significant parts of this article rests on //'-like norms. Largely, 
this choice is motivated by the numerical scheme, namely, an //' -conforming FEM. 



1.1 Non-standard FEM 

The limitations of the classical FEM mentioned above in (HJl and (|ll| have sparked a 
significant amount of research in the past decades to overcome or at least mitigate 
them. This research focuses on two techniques that are often considered in tandem: 
firstly, the underlying approximation by classical piecewise polynomials is replaced 
with special, problem-adapted functions such as systems of plane waves; secondly, 
the numerical scheme is based on a variational formulation different from the classi- 
cal //'-conforming Galerkin approach. Before discussing these ideas in more detail, 
we point the reader to the interesting work |1T2], which shows for a model situation 
on regular, infinite grids in 2D that no 9-point stencil (i.e., a numerical method based 
on connecting the value at a node with those of the 8 nearest neighbors) generates a 
completely pollution-free method; the ID situation is special and discussed briefly 
in Section Q 

Work that is based on a new or modified variational formulation but rests on 
the approximation properties of piecewise polynomials includes the Galerkin Least 
Squares Method II39II40I . the methods of ID, and Discontinuous Galerkin Methods 
( II33H35I and references there). Several methods have been proposed that are based 
on the approximation properties of special, problem-adapted systems of functions 
such as systems of plane waves. In an //'-conforming Galerkin setting, this idea has 
been pursued in the Partition of Unity Method/Generalized FEM by several authors, 
e.g., |l5l|45l|50l|5l]|56l|60l|68l|69l|8ll- A variety of other methods that are based on 
problem-adapted ansatz functions leave the //'-conforming Galerkin setting and en- 
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force the jump across element boundaries in a weak sense. This can be done by least 
squares techniques ( II 1 4115 3 165ll70ll80l and references there), by Lagrange multiplier 
techniques as in the Discontinuous Enrichment Method r31"32!,l821 or by Discon- 
tinuous Galerkin (DG) type methods, [19-21 ,36, 42, 43 , 46, 55,63, 64|; in these last 
references, we have included the work on the Ultra Weak Variational Formulation 
(UWVF) since it can be understood as a special DG method as discussed in II19II36L 



1.2 Scope of the article 

The present article focuses on the stability properties of numerical methods for 
Helmholtz problems and exemplarily discusses three different approaches in more 
detail for their differences in techniques. The first approach, studied in SectionlH is 
that of the classical //^ -based Galerkin method for Helmholtz problems. The setting 
is that of a Garding inequality so that stability of a numerical method can be inferred 
from the stability of the continuous problem by perturbation arguments. This moti- 
vates us to study for problem (|9]l, which will serve as our model Helmholtz problem 
in this article, the stability properties of the continuous problem in Section |2] In 
order to make the perturbation argument explicit in the wavenumber k, a detailed, k- 
explicit regularity analysis for Helmholtz problems is necessary. This is worked out 
in Section |3]for our model problem (|9) posed on polygonal domains. These results 
generalize a similar regularity theory for convex polygons or domains with analytic 
boundary of II61II62I . Structurally similar results have been obtained in connection 
with boundary integral formulations in 05411591 . 

We discuss in Sections 16.21 and 16.31 somewhat briefly a second and a third ap- 
proach to stability of numerical schemes. In contrast to the setting discussed above, 
where stability is only ensured asymptotically for sufficiently fine discretizations, 
these methods are stable by construction and can even feature quasioptimality in 
appropriate residual norms. We point out, however, that relating this residual norm 
to a more standard norm such as the L^-norm for the error is a non-tiivial task. Our 
presentation for these methods will follow the works |fT9ll36ll43ll65l . 

Many aspects of discretizations for Helmholtz problems are not addressed in this 
article. For recent developments in boundary element techniques for this problem 
class, we refer to the survey article [|22l . The model problem ^ discussed here 
involves the rather simple boundary condition (|9bi . which can be understood as 
an approximation to a Dirichlet-to-Neumann operator that provides a coupling to 
a homogeneous Helmholtz equation in an exterior domain together with appropri- 
ate radiation conditions at infinity. A variety of techniques for such problems are 
discussed in 1371 . Further methods include FEM-BEM coupling, the PML due to 
Berenger (see III7II24I and references therein), infinite elements ll26l . and methods 
based on the pole condition, ll44l . Another topic not addressed here is the solution 
of the arising linear system; we refer the reader to [I28JI301 for a discussion of the 
state of the art. Fuither works with survey character includes Il29ll47ll48l[83l . 
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We employ standard notation for Sobolev spaces, |[T][T8]|67]|221 . For domains co and 
ky^Qwe denote 

\\4lt.--=kH4h^o>) + \\'^4l2(o,y (7) 

This norm is equivalent to the standard H ' -norm. The presence of the weight k in the 
L'^-part leads to a balance between the H^-seminorm and the L^-norm for functions 
with the expected oscillatory behavior such as plane waves e''^^'" (with d being a unit 
vector). Additionally, the bilinear form B considered below is bounded uniformly in 
k with respect to this (A:-dependent) norm. 

Throughout this work, a standing assumption will be 

1^1 > ^0 > 0; (8) 

our frequently used phrase "independent of k" will still implicitly assume ([8]l. We 
denote by C a generic constant. If not stated otherwise, C will be independent of 
the wavenumber k but may depend on k(). For smooth functions u defined on a d- 

dimensionalmanifold, we employ the notation I V"m(x)|^ ;= ^ 

aeNn: \a\=n 



1.4 A model problem 

In order to fix ideas, we will use the following, specific model problem: For a 
bounded Lipschitz domain Q. cW', d e {2,3}, we study for e K, > k^, the 
boundary value problem 

- Aii-k^u ^ fxwQ., (9a) 
d„u + iku ^ gon do. . (9b) 

Henceforth, to simplify the notation, we assume A: > A:o > but point out that the 
choice of the sign of k is not essential. The weak formulation for (|9]l is; 

FindMe//i(i2) s.t. B{u,v) = l{v) VvG//'(i2), (10) 

where, for / G L^(i2) and g e L^{dQ.), B and I are given by 

Biu,v):= JjVu-Vv-k\v)+ik J^^uv, l{v):={f,v)ii(a) + {8^^)L^dny (H) 

As usual, if / e {H^ {Q.))' and g € H-^/-{dn), then the -inner products ( • , • )i^2 /Q-t 
and (•, OiSf^ij) are understood as duality pairings. The multiplicative trace inequal- 
ity proves continuity of B; in fact, there exists Cb > independent of k such that 
(see, e.g., IMl Cor 3.2] for details) 
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2 Stability of the continuous problem 

Helmholtz problems can often be cast in the form "coercive + compact perturbation" 
where the compact perturbation is A:-dependent. In other words, a Garding inequality 
is satisfied. For example, the sesquilinear form B of ( fTTT i is of this form since 



and the embedding H^{£2) C L^{£2) is compact by Rellich's theorem. Classical 
Fredholm theory (the "Fredholm alternative") then yields unique solvability of dTol i 
for all / e (//' (Q))' and g e H-^l^{dQ.), if one can show uniqueness. Uniqueness 
in turn is often obtained by exploiting analyticity of the solutions of homogeneous 
Helmholtz equation, or, more generally, the unique continuation principle for elliptic 
problems, (see, e.g., Il52l Chap. 4.3]): 

Example 2.1 (Uniqueness for (|9)). Let / = and g = in Then, any solu- 
tion u e H^{Q) of (|9]l satisfies u\gQ = since = ImB{u,u) = ^Il"ll^2(5i2) 
Lemma |2!2] |. Hence, the trivial extension m to satisfies u £ (H^^)- The obser- 
vations B{u,v) = for all v E H^{Q) and u\gQ =0 show 



Hence, m is a solution of the homogeneous Helmholtz equation and u vanishes on 
\ f2. Analyticity of u (or, more generally, the unique continuation principle pre- 
sented in 1521 Chap. 4.3]) then implies that u = 0, which in turn yields m = 0. ■ 

The arguments based on the Fredholm alternative do not give any indication of how 
the solution operator depends on the wavenumber k. Yet, it is clearly of interest 
to know how k enters bounds for the solution operator It turns out that both the 
geometry and the type of boundary conditions strongly affect these bounds. For ex- 
ample, for an exterior Dirichlet problem, |116 1 exhibits a geometry and a sequence of 
wavenumber (A;„)„gp} tending to infinity such that the norm of the solution operator 
for these wavenumbers is bounded from below by an exponentially growing term 
Ce''''" for some C, b > 0. These geometries feature so-called "trapping" or near- 
trapping and are not convex. For convex or at least star-shaped geometries, the k- 
dependence is much better behaved. An important ingredient of the analysis on such 
geometries are special test functions in the variational formulation. For example, as- 
suming in the the model problem (fTOl i that Q is star-shaped with respect to the ori- 
gin (and has a smooth boundary), one may take as the test function v(x) = x ■ Vm(x), 
where u is the exact solution. An integration by parts (more generally, the so-called 
"Rellich identities" 1671 p. 261] or an identity due to Pohozaev, fTTl ) then leads to 
the following estimate for the model problem (ITOl l: 



RsB{u,u) +2k^{u,u)^HQ-j ^ hWik.n 



(13) 
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||M||ii',r2 <C \\f\\L2(Q) + \\8\\L^(dn) ; 



(14) 



this was shown in 1561 Prop. 8.1.4] (for d = 2) and subsequently by fZ5\ for d = 3. 
Uniform in k bounds were established in BTl for star-shaped domains and certain 
boundary conditions of mixed type by related techniques. The same test function 
was also crucial for a boundary integral setting in ll23l . A refined version of this 
test function that goes back to Morawetz and Ludwig, ll66l was used recently in a 
boundary integral equations context (still for star-shaped domains), ll78l . 

While (fl4] i does not make minimal assumptions on the regularity of / and g, the 
estimate (fT4l i can be used to show that (for star-shaped domains) the sesquilinear 
form B of dTol i satisfies an inf-sup condition with inf-sup constant 7 = 0{k^ ' ) — this 
can be shown using the arguments presented in the proof Theorem |2.5l 

An important ingredient of the regularity and stability theory will be the concept 
of polynomial well-posedness by which we mean polynomial-in-fc-bounds for the 
norm of the solution operator The model problem (|9) on star- shaped domains with 
the a priori bound (fl4l i is an example. The following Section lTTI shows polynomial 
well-posedness for the model problem (|9]l on general Lipschitz domains (Thm. |24b . 
It is thus not the geometry but the type of boundary conditions in our model problem 
(|9]l, namely, Robin boundary conditions that makes it polynomially well-posed. In 
contrast, the Dirichlet boundary conditions in conjunction with the lack of star- 
shapedness in the examples given in lfT6l make these problem not polynomially 
well-posed. 



2.1 Polynomial well-posedness for the model problem ® 

Lemma 2.2. Let Q, C be a bounded Lipschitz domain. Let u G Ll^ (Q) be a weak 
solution of ^ with f ^ and g £ L^{d £2). Then \\u\\[2(^gQ^ ^ k^^\\8\\L^{dn)- 

Proof. Selecting v = m in the weak formulation (fTOl i and considering the imaginary 
part yields 



This concludes the argument. □ 

Next we use results on layer potentials for the Helmholtz equation from ll59l to 
prove the following lemma: 

Lemma 2.3. Let £2 C be a bounded Lipschitz domain, m G //' (i2) solve (|2| with 
/ = 0. Assume mI^^j G L^idQ.) and d„u G L^{dQ). Then there exists C > indepen- 
dent ofk and u such that 





\WnU\\i2(s^) . 
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Proof. With the single layer and double layer potentials V^. and we have the 
representation formula u ~ VkdnU~KkU- From |f59l Lemmata 2.1, 2.2, Theorems 4.1, 
4.2] we obtain 



Thus, 



MiHn) < Ck (||M||i2(5^) + \\d„u\\H-i^dn) 
Next, using v = u in the weak formulation (fTOl l yields 



IVmI 



<C 



^^ll"llL2(i2) + ll^"«llz.2(ai2)ll«IL2(ai2) 



and therefore 



IVmII 



r,2||,,i|2 



-k^UM 



r,-4| 



H-^{dn) ' 



which concludes the proof. □ 

Theorem 2.4. Let £2 C M'', d e {2,3} be a bounded Lipschitz domain. Then there 
exists C > (independent of k) such that for f £L?'{£2) and g G L?{dQ.) the solution 
ueH^Q) o/lO satisfies 



Ml.k.Q < c 



'k%\\LHdn)+k'^\f\\LHn) 



Proof. We first transform the problem to one with homogeneous right-hand side / 
in the standard way. A particular solution of the equation (|9a] i is given by the Newton 
potential mq : = Gj^ here, is a Green's function for the Helmholtz equation and 
we tacitly extend / by zero outside Q. Then uq G Hf^^^lW') and by the analysis of 
the Newton potential given in 1611 Lemma 3.5] we have 



k '!l«o||//2(r2) + ll«o||//l(i^) +^||Mo|lL2(i2) < C||/||^2(^2) 

The difference u: — u — uq then satisfies 

-Au-k^u = Q mQ, 
d„u + iku ~ g — (<3„Mo + ikuo) =: g. 

We have with the multiplicative trace inequaUty 

WsWLHdn) < C lkllL2(ar2) + ll"o||]y'?(i2)ll"oll//Hi2)+^ll"ollHi^(i2)l 
<c[\\gh2(dn)+k'^^\\f\\LHn)_ ■ 
To get bounds on u, we employ Lemma |2T2l and ( [TtI i to conclude 



(15) 



(16a) 
(16b) 



I«0|li2(i2) 



(17) 
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Mmdn) < Ck-'\\g\\i2^sn) < C [k-'\\g\\LHda) +k-''^\\f\\LHn]] , (18) 
WdrMtHda) < c\\\g\\LHda)+k\\'A\LHda)\ <c\\\g\\ii^;)a)+k^'^\\f\\LHa)\- (19) 
Lemma lZll and the generous estimate \\dnu\\}j-i(dn) — C|l^n"llL2(3i2) produce 

(20) 

Combining ( fTSl ). ( l20b finishes the argument. □ 

The a priori estimate of Theorem 12.41 does not make minimal assumptions on the 
regularity of / and g. However, it can be used to obtain estimates on the inf-sup and 
hence a priori bounds for / e (//' (i2))' and g e H^^I^{dQ.) as we now show: 

Theorem 2.5. Let Q. c W', d ^ {2^3} be a bounded Lipschitz domain. Then there 
exists C > (independent of k) such that the sesquilinear fonn B of Ml^ satisfies 



inf 



sup 



ReB(M,v) 



> ck-y\ 



(21) 



Furthermore, for every f <E (//'(i2))' and g ^ H '/^((9i2) the problem ( TTOl ) is 
uniquely solvable, and its solution u £ H^{£2) satisfies the a priori bound 



\u\\i,k.n < Ck 



7/2 



\\f\\{Hi{n)y + \\8\\H-'P-{dn] 



(22) 



If Q is convex or if Q is star-shaped and has a smooth boundary, then the following, 
sharper estimate holds: 



ReB{u,v) 
mf sup — — 



>Ck 



-1 



(23) 



Proof. The proof relies on standard arguments for sesquilinear forms satisfying a 
Garding inequality. For simplicity of notation, we write || • || i.i for || • || i,jt,i2. 
Given u G H^{Q) we define zG H^{Q) as the solution of 

2^^(-,m)l2(i2) ='B(-,z). 

Theorem l2.4l implies \\z\\i,k < C^^/^ j|M||^2(^), and v — u + z satisfies 



ReB{u,v) ^RcB{u,u)+ReB{u,z) 
Thus, 



l.k ' 



2k^\\u\\l2^^^+ReB{u,z) = \\u\\l,. 



ReB{u,v) = Mil,, 

hhk = ||«+z||l,A-< ||M||l,i+||z||ll-< ||Mi|l,i + C^''/^||MilL2(I2) <Cfc'''^l|M|ll,^ 
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Therefore, 



ReB(M,v) = |l«|l2,>||„||i,,C^-7/2||^,||j^^^ 



which concludes the proof of (ISTT l. Example l2. 1 I provides unique solvability for (|9]l 
so that dJTJi gives the a priori estimate ( |22] |. Finally, ( |23] | is shown by the same 
arguments using (fl4] i. □ 



3 ^-explicit regularity theory 
3.1 Regularity by decomposition 

Since the Sobolev regularity of elliptic problems is determined by the leading order 
terms of the differential equation and the boundary conditions, the Sobolev regu- 
larity properties of our model problem (0 are the same as those for the Laplacian. 
However, regularity results that are explicit in the wavenumber k are clearly of in- 
terest; for example, we will use them in Section |4j2]below to quantify how fine the 
discretization has to be (relative to k) so that the FEM is stable and quasi-optimal. 

The A:-explicit regularity theory developed in [16 111621 (and, similarly, for integral 
equations in 15411591 ) takes the form of an additive splitting of the solution into 
a part with finite regularity but ^-independent bounds and a part that is analytic 
and for which fc-explicit bounds for all derivatives are available. Below, we will 
present a similar regularity theory for the model problem (0) for polygonal X2 C M^, 
thereby extending the results of ll62l . which restricted its analysis of polygons to the 
convex case. In order to motivate the ensuing developments, we quote from f6l\ a 
result that shows in a simple setting the key features of our ^-explicit "regularity by 
decomposition": 

Lemma 3.1 ( tSB Lemma 3.5J). Let Br{0) dW'.de {1,2,3} be the ball of radius 
R centered at the origin. Then, there exist C, 7 > such that for all k (with k > ko) 
the following is true: For all f € L^(M'') with supp/ C Br{Q) the solution u of 



Au-k^u^f inW' , 



subject to the Sommerfeld radiation condition 




for \x\ -> 



00 



has the following regularity properties: 
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^II"//2|Ii,M2k(0) + II"//2|I//2(B2r(0)) - ^II/IIl2(B;;(0)), 

I|V"«-c/IIl2(b2,(0)) < Cy" max{«,n""'||/IL2(B,(0)) V« e No- 

A few comments concerning Lemma ITTI are in order For general / e L^(B^(0)), 
one cannot expect better regularity than //^-regularity for the solution u and, in- 
deed, both regularity results (|i} and dn} assert this. The estimate ^ is sharp in 
its dependence on k as the following simple example shows: For the fundamen- 
tal solution Gk (with singularity at the origin) and a cut-off function x G Cq{W^) 
with supp;|^ C B2r{Q) and = 1 on Br{Q), the functions u := (1 — x)Gk and 
/ -Au-k^u satisfy ||m||h2(b,^(o)) = 0{k^) and ||/!li2(B^(o)) = 0{k). Compared 
to dill, the regularity assertion dn) is finer in that its //^-part u^i has a better k- 
dependence. The ^-dependence of the analytic part m.c/ is not improved (indeed, 
II"j^II//2(b2r(0)) ^ C^II/IIl2(Br(0)))' but the analyticity of is a feature that higher 
order methods can exploit. 

The decomposition in dUll of Lemma [TT] is obtained by a decomposition of the 
datum / using low pass and high pass filters, i.e., / = L^i^f + H^i^f, where the low 
pass filter L^j^ cuts off frequencies beyond r\k (here, T] > 1) and //^<. eliminates 
the frequencies small than 7]k. Similar frequency filters will be important tools in 
our analysis below as well (see Sec. 13.3.1b . The regularity properties stated in dH) 
then follow from this decomposition and the explicit solution formula m = G^*/ 
(see IMI Lemma 3.5] for details). 

Lemma lTTI serves as a prototype for "regularity theory by decomposition". Sim- 
ilar decompositions have been developed recently for several Helmholtz problems 
in ll62l and ll54ll59l . Although they vary in their details, these decomposition are 
structurally similar in that they have the form of an additive splitting into a part with 
finite regularity with ^-independent bounds and an analytic part with fc-dependent 
bounds. The basic ingredients of these decomposition results are (a) (piecewise) 
analyticity of the geometry (or, more generally, the data of the problem) and (b) a 
priori bounds for solution operator The latter appear only in the estimate for the 
analytic part of the decomposition, and the most interesting case is that of polyno- 
mially well-posed problems. We illustrate the construction of the decomposition for 
the model problem dU in polygonal domains f2 C K^. This result is an extension 
to general polygons of the results ll62l . which restricted its attention to the case of 
convex polygons. We emphasize that the choice of the boundary conditions d9b] i is 
not essential for the form of the decomposition and other boundary conditions could 
be treated using similar techniques. 



3.2 Setting and main result 

Let X2 C K'^ be a bounded, polygonal Lipschitz domain with vertices Ay, y = 1 , . . . , 7, 
and interior angles Oj, 7 = 1, ... ,7. We will require the countably normed spaces 
introduced in ||8.57J . These space are designed to capture important features of solu- 



On stability of discretizations of the Helmholtz equation (extended version) 



13 



tions of elliptic partial differential equations posed on polygons, namely, analyticity 
of the solution and the singular behavior at the vertices. Their characterization in 
terms of these countably normed spaces also permits proving exponential conver- 
gence of piecewise polynomial approximation on appropriately graded meshes. 

These countably normed spaces are defined with the aid of weight functions 
^p-^ 1^ that we now define. For j3 e [0, 1), « eNo,k> 0, and j3 G [0, 1)'', we set 

u+j3 

Irl I 

^«./3.A:W = min<( 1, 



min 



^„tAW = ri^«,P„^(--^;)- (24) 

7=1 

1 /2 

Finally, we denote by Hp„ {dQ) the space of functions whose restrictions of the 
edges of dQ are in //'/^. 

We furthermore introduce the constant Csoi{k) as a suitable norm of the solution 
operator for (|9). That is, Csoi{k) is such that for all f <E L?{Q.), g <E L^{dQ) and 
corresponding solution u of (|9]) there holds 



u\\i.Ln < Qoi{k) ||/|lL2(i2) + WshndQ) 



(25) 



We recall that Theorem 12.41 gives Csoi{k) = 0{k^/^) for general polygons and 
Csoi{k) = 0{1) by ll56l Prop. 8.1.4] for convex polygons. Our motivation for using 
the notation Csoi{k) is emphasize in the following theorem how a priori estimates 
for Helmholtz problems affect the decomposition result: 

Theorem 3.2. Let £2 C be a polygon with vertices Aj, j — 1,...,J. Then there 
exist constants C, y > 0, P [0,1)"' independent of k > ko such that for every 
f E L^{Q) and g € Hpl^idCl) the solution u of ^ can be written as u~ Uj^i +Uj^ 
with 

^||M//2!|L/:,i2 + ||M//2|1//2(^2) < CC/,g 

ll«^/ll//i(^2) - iC.U,l{k) + l)Cf^g 

klW^hnn) < {Csoi{k)+k)Cf,j, 
Il^„.^,^^"^'«-^^lli^(i2) ^ C{CUk) + l)Ar-'max{«,Ar}"+V'C/,, V« S No 

with Cf,g ||/||l2(^) + ||g||„i/2,;,„> andCsoi{k) introduced in ^25^ . 



Proof. The proof is relegated to Section 13.41 We mention that the A:-dependence of 
our bounds on ||M.c/||i2(^) is likely to be suboptimal due to our method of proof. □ 



Theorem 13.21 mav be viewed as the analog of Lemma I3TTI dnji; we conclude this 
section with the analog of Lemma ITTl 
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Corollary 3.3. Assume the hypotheses of Theorem 13.21 Then there exist constants 
C>Q,~^ e [0,iy independent of k such that for all f e L^{Q), g eHp(^{dQ) the 
solution u o/dP]) satisfies \\u\\ ij^q ^ CCsoi{k) \\f\\]2 



as well as 



\\%,^,^MLHn)<Ck{Qo,{k) + l) \\f\\LHn) + M„;j2^,^^ 



Proof. The estimate for ||m|| i q expresses dZSl l. The estimate for the second deriva- 
tives of u follows from Theorem 13.21 since u = Ujji + Uj;/. □ 



3.3 Auxiliary results 

Just as in the proof of Lemma 13.11 an important ingredient of the proof of Theo- 
rem [3]2] are high and low pass filters. The underlying reason is that the Helmholtz 
operator —A — k^ acts very differently on low frequency and high frequency func- 
tions. Here, the dividing line between high and low frequencies is at 0{k). For 
this reason, appropriate high and low pass filters are defined and analyzed in Sec- 
tion 13.3.11 Furthermore, when applied to high frequency functions the Helmholtz 
operator behaves similarly to the Laplacian —4 or the modified Helmholtz operator 
—A +k^. This latter operator, being positive definite, is easier to analyze and yet 
provides insight into the behavior of the Helmholtz operator restricted to high fre- 
quency functions. The modified Helmholtz operator will therefore be a tool for the 
proof of Theorem 122] and is thus analyzed in Section l3.3.3l 



3.3.1 High and low pass filters, auxiUary results 

For the polygonal domain i2 C we introduce for 7] > 1 the following two low 
and high pass filters in terms of the Fourier transform JF: 

1. The low and high pass filters LQ^r)f ■■ L^i^) ^ L^{^) and Hq,j^ : L^(Q) 
L^{Q) are defined by 

Ln.^f^{-^-'XB,,(o)^iEnf))\n, Hn.^f ^ {-^-'xM^XB^^ioy^iEnfm; 

here, Bj^kiO) is the ball of radius rjk with center 0, the characteristic function of 
a set A is Xa^ and Eq denotes the Stein extension operator of ||79] Chap. VI]. 

2. Analogously, we define L^n^r^f : L^{dQ) L^{dQ) and H^q j^ : L^{dQ) 

{dQ ) in an edgewise fashion. Specifically, identifying an edge e of Q with 
an interval and letting Eg be the Stein extension operator for the interval e C R 
to the real line R, we can define with the univariate Fourier transformation 
the operators Li,,n and Hf,n by 
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LeM={'^''XB,M'^{E,g))\e, He.^g={^-'xm\B,M-^{EJ))\e; 

the operators and//^^ ,j are then defined edgewise by {L^Q.qg)\e~Le^r\g 

and {Hdn^r\g)\e = He,^g for all edges ecdQ. 

These operators provide stable decompositions of {Q ) and L?{dQ ). For example, 
one has Lq +Hq j^ = Id on L^{Q) and the bounds 

where C > depends solely on Q (via the Stein extension operator £^2)- The op- 
erators Hq ji and HgQ jj have furthermore approximation properties if the function 
they are applied to has some Sobolev regularity. We illustrate this for H^q : 

Lemma 3.4. Let Q <Z M? be a polygon. Then there exists C > independent of k 

1 /2 

and 77 > 1 such that for all g G Hp„ {d£2 ) 



H''t{dn) 



Proof. We only show the estimate for \\Hdn,r]g\\ii(gQy We consider first the case of 
an interval/ C M. We define Hi,-^g by Hi^^^g ^ ^''^XM.\Br,,,(Q)-^Eig, where Zm\b^j.(0) 
is the characteristic function for R\ (—17 A:, r\k) and £/ is the Stein extension operator 
for the interval /. Since, by Parseval, is an isometry on (R) we have 

|2 



\\Hi,ng\\i2(j) < \\Hi.r^g\\i2^^) = i „ _ \^E,g\'d^ 



The last integral can be bounded by Clj/i/gll^i/j^jj)- The stability properties of the 
extension operator £/ then imply furthermore ||£'/g||^i/2(]u) < C||g||^i/2(^) . In total. 



\H,. 



where, in the last estimate, the constant C depends additionally on kQ. From this 
estimate, we obtain the desired bound for |l//ai2,r]§llL2(ai2) identifying each edge 
of Q with an interval. □ 



3.3.2 Corner singularities 



We recall the following result harking back to the work by Kondratiev and Grisvard: 
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Lemma 3.5. Let £1 c be a polygon with vertices Aj, j ~ 1, ... ,7, and interior 
angles COj, j = 1 , ... ,7. Define for each vertex Aj the singularity function Sj by 

5,(o,(p,,) = r;/"^cos(^|-(p,-), (26) 

where {rj,(pj) are polar coordinates centered at the vertex A j such that the edges of 
£2 meeting at Aj correspond to (pj ~ and (pj = COj. Then every solution u of 

—Au^f in Q, d„u = g ondQ, 

can be written as u ^ uq + Lj=i '^f if:8)^j '^^th the a priori bounds 

\Wo\\H2in) + t\^f(f'8)\<C 
j=i 

The a^ are linear functionals, and a^ = for convex corners A j (i.e., ifcOj < n). 
Proof. This classical result is comprehensively treated in Il38l . □ 



\f\y 



(27) 



3.3.3 The modified Helmlioltz equation 

We consider the modified Helmholtz equation in both a bounded domain with Robin 
boundary conditions and in the full space M^. The corresponding solution operators 
will be denoted 5^ and S^2- 

1. The operators'^ : L^{Q) x Hpw{dQ) H^{Q) is the solution operator for 

— Au + k^u = f ini2, d„u + iku = g ondQ. (28) 

2. The operator ■ L^{R^) (R^) is the solution operator for 

-Au + k^u^f inM^. (29) 

Lemma 3.6 (properties of 5^). Let Q <zM? be a polygon and f £ L^^Q), g & 
Hpw {d£2). Then the solution u := SQ{f,g) satisfies 

Mixn < k-'^'\\g\\L2idn)+k-'\\f\\LHny (30) 

Furthermore, there exists C > independent ofk and the data f, g, and there exists 
a decomposition u = m^2 +Y4=iCi^{f,g)Sifor some linear functionals aj^ with 

J 

il«//2|l//2(i^) + El«^(/'^)l 

1=1 



ll/!lL2(i^) + llgll^V2(5^)+^''''kllL20i2) 



(31) 
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Proof. The estimate ( [30l l for 1 1 m 1 1 1 follows by Lax-Milgram - see 1621 Lemma 4.6] 
for details. Since u satisfies 

—Au = f — k^u in f2 , d„u ^ g — iku on dQ , 

the standard regularity theory for the Laplacian (see Lemma 13.51 ) permits us to de- 
compose M = Uff2 + Y/i=i if ^ k^u,g — iku)Si. The continuity of the linear func- 
tionals af reads 

J 

Y,\afif~k\,g~iku)\<C 

1=1 

Since {f,g) ^ is linear, the map {f,g) a+{f,g) -.^ af {f ~ k^u,g -iku) 

is linear, and ^0^, (|27]i give the desired estimates for m^2 and (/j,?)- n 

Lemma 3.7 (properties of 5j|^2)- There exists C > such that for every 77 > 1 and 
every f £ L^{M?) whose Fourier transform ^/ satisfies suppj^/ C \ BrikiO), 
the solution u = ^^^2/ o/(|29l) satisfies 

Proof. The result follows from Parseval's theorem and the weak formulation for u 
as follows (we abbreviate the Fourier transforms by / ^f and u — ^u): 

II"IIi,A-.M2 (/:«)l2(M2) = (/7")l2(M2) 

= J , (l^l^+^^)-M/P^^H«lll...«2 < , , ^ J l/llL2(R2)|k^||i,,,K2, 

where, in the penultimate step, we used the support properties of /. Appealing again 
to Parseval, we get the desired claim for j|M|| 1 r2. The estimate for ||m||//2(jj2) now 
follows from f £ (R^) and the standard interior regularity for the Laplacian. □ 



ll/-^'«ilL2(I2) + k-ifc«ll^l/2(5^) • 



3.4 Proof of Theorem 13.21 

We denote by S : {f,g) S{f,g) the solution operator for (|9}. Concerning some of 
its properties, we have the following result taken essentially from ||62l Lemma 4.13]: 

Lemma 3.8 (analytic regularity of S{f,g)). Let Q be a polygon. Let f be analytic 
on £2 and g £ L^ [dO, ) be piecewise analytic and satisfy for some constants Cf, Cg, 



18 S. Esterhazy and J.M. Melenk 
7/. 7g > 

I|V'7IIl2(^) <C///max{«,n« Vn G No (32a) 

l|V'kllL2w <Q^'max{«,fc}" VneNo Ve G <f , (32b) 

where S denotes the set of edges of Q. and Vj- tangential differentiation. Then there 
exist /3 G [0, 1 f depending only on Q) and constants C, 7 > f depending only on 
Q, Yf, Yg, ko) such that the following is true with the constant Csoi{k) of ( 1251 ); 

\\4i.k.n <Qoi{k){Cf + Cg) (33) 
^ CQ,i{k)k-\Cf + Cg)fmax{n,kr+^ \fn e Nq. (34) 

Proof. The estimate ( |33] | is simply a restatement of the definition of Csoi{k). The 
estimate ( l34l i will follow from ll57l Prop. 5.4.5]. To simplify the presentation, we 
assume by linearity that g vanishes on all edges of £2 with the exception of one edge 
r. Furthermore, we restrict our attention to the vicinity of one vertex, which we take 
to be the origin; we assume F C (0,°°) x {0}, and that near the origin, Q is above 
(0,°o) X {0},i.e., {(rcos(p,rsin(p): 0<r<p,0<(j!)<aj}ci2 for some p,(0>Q. 
Upon setting e := we note that u solves 

—£^Au — u~e^f oni2, £^d„u = £{£g — iu) ondil. 

On the edge F, the function g is the restriction of Gix){x,y) := g{x)e^^'^^ to F. The 
assumptions on / and g then imply that 1571 Prop. 5.4.5] is applicable with the 
following choice of constants appearing in | |57l Prop. 5.4.5]: 

Cf = £^Cf, ~ ee^/^Cg, Cc, = e, Cj, = 0, = 1, 
7/ = 0(1), 7g. =0(1), 7g;=0(1), 7i = 0, 7, = 0, 

resulting in the existence of constants C,K>0 and j8 e [0, 1)'' with 

^ Cr'+2max{« + 2,n"+' (^-^C/ + ^-'||«|| 1.,,^ 

for all n G Nq. We conclude the argument by inserting (|3?t and estimating gener- 
ously k-^Cf + k- ^I^Cg <c{Cf + . 

We remark that this last generous estimate comes from the precise form of our 
stability assumption ( |25] l. Its form ( |25] l is motivated by the estimates available for 
the star-shaped case, but could clearly be replaced with other assumptions. □ 

Corollary 3.9 (analytic regularity of S{Lq j^f,LsQ ng)). Let Q be a polygon and 

TJ > 1. Then there exist j5 G [0, 1 )'' f depending only on Q ) and C, 7 > f depending 
only on £2, ko, and rj > I) such that for every / G L^{£2) and g G L?{dQ), the 
junction u = SiLn^y^f ,L3Q r\g) satisfies with Cf^g := ||/||i2(^) + 11^11^2(5^2) 
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\\uh.k.a<CC,oi{k)Cf.s (35) 
||f„^,V"+2„||^2(^) <CC,w(^)^"V'max{«,^}"+2C/,, V« e Nq. (36) 

Proof. The definitions of LQnf and L^q imply with Parseval 

l|V"Lr2,r,/|lz,2(^2) < C(77^r!l/IL2(^) V« e No, 
l|V"rLar2,,g|L2(,^) < C(77fc)«||g||^2(,^) Vn e No, 

where again Vj- is the (edgewise) tangential gradient. The desired estimates now 
follow from Lemma[378] □ 



Key to the proof of Theorem l3.2l is the following contraction result: 
Lemma 3.10 (contraction lemma). Let £2 ciR^ be a polygon. Fix q e (0, 1). Then 
one can find j3 G [0, !)■' (depending solely on Q ) and constants C, y > indepen- 
dent ofk such that for every / G L^(f2) and every g £ Hpw{dQ), the solution u of 
(O can be decomposed as u = Ufji + \ ai{f,g)Si + + r, where Uffi € H^{Q), 
the at are linear functionals, anduj^ G C°°(X2). These functions satisfy 

J 



k\\uHA\\,k.n + Il«//2|l//2(i2) + L < C Il/|lz,2(i2) - 

1=1 

\\u.A\,k..n<CCsoi{k) ||/||l2(^2) + 11^11^2(5^2) ' 

< CC„„W^- Vmax{«,n"+' [||/|Il2( 
for all n e No- Finally, the remainder r satisfies 

—Ar — k^r = f onQ, d„r + ikr — g 
for some f e L^{Q) andge Hl,i^{dQ) with 



Proof We start by decomposing {f,g) = [L^.r^f^LdQ-^g) + {Ha.j^fMgn.rig) with 
a parameter 77 > 1 that will be selected below. We set 



Mi/ S{LQ„f,LgQ . 



where we tacitly extended H^ -qf (which is only defined on Q) by zero outside 
12. Then Mj,/ satisfies the desired estimates by Corollary 13.91 For u\ we have by 
Lemma [3^ and the stability \\Ha_r\f\\i2(^Q-^ < C||/||^2(^) (we note that C > is 
independent of k and rj ) the a priori estimates 



II 1,,,M2 <C^-Hl + 772)-l/2||//^ ,,/||^,^^j<C^-l(l+^)-l 11/11 

|mi||//2(r2.| < C\\Hn,iTf\\i2i^n) < CWfW^^i^^y 
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The trace and the multiplicative trace inequalities imply for^i := ^„mi +ikui: 

k'^Hl + rin8i\\LH3n) + \\8i\\^in^^^^ 
For g2 := H;jQ j^g — we then get from Lemma 13741 and the triangle inequality 



Lemma [T6] vields for U2 '■= {Q,g2), 

\\u2\\i.k.n<Ck-'/^\\g2h2idQ)<Ck-\l+il)-''^ 



l/llL2(i2) + l|g||„l/2(^^) 



and furthermore we can write U2 — u^i +L/=i ^i{^Tg2)Si, with 



(0,^2)1 <C 



l/llL2(r2) + kll^V20i2) 



We then define ai{f,g) := 0,^(0,^2) and note that {f,g) n- ai{f,g) is Unearby lin- 
earity of the maps aj' and {f,g) 1— >■ g2- The above shows that u^i and the a,- satisfy 
the required estimates. Finally, the function u:^u — {uj^/ + ui+ U2) satisfies 



-Au — k^u = 2k^(ui + U2) —'■ /, dnU + iku = =: , 



together with 



\\f\\L2(Q)<2k^\\uih2,n) + h2\\L2in)) <C(l + 7? 



-1/2 



l/ll 



L2{n) + \\8\ 



Hence, selecting 77 > 1 sufficiently large so that for the chosen <7 e (0, 1) we have 
C{1 +rj)^^^~ < q allows us to conclude the proof. □ 



Proof of Theorem 13.21 The contraction property of Lemma 13.101 can be restated 
as S{f,g) = Uff2 +T/i=iai{f,g)Si + u_^ + S{f,g), where, for a chosen q G (0,1), 



we have ||/||i2(^) + \\g\\ 



< 



together with ap- 



11/11^2(12)- 

propriate estimates for m^2, ai{f,g), and Uj^. This consideration can be repeated 
for S{f,g). We conclude that a geometric series argument can be employed to write 
u = S{f,g) = Uff2 + Y/i=i^i{f,8)Si + where m^2 G H^{Q), m_0/ e C°°(i2), and 
the coefficients a, are in fact linear functionals of the data {f,g)- Furthermore, we 
have with the abbreviation C/,g ;= 11/11^2(^2) + kll^^i/ij^^) 
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\Wv/\\ i.k.n < CCf,g 

J 

i=l 

Finally, Lemma |3 . 11 1 below allows us to absorb the contribution Lf=i in 
the analytic part by setting := + Y/i=i '^i{f^g)Si- In view of j3,- < 1, we have 
2 — j3,- > 1 and arrive at 

hAuHn) < C{C,oi{k) + l)C/,g, ^ll"i'/llL2(i2) < CC/,g(C,„,(^) +^), 
^ '^'^/.s +^"'] max{«,fc}"+2 V« e No, 

which concludes the argument. □ 

Lemma 3.11. Let j3; e [0, 1) satisfy )3,- > 1 - ^. Then, for some C, 7> independent 
ofk, the singularity functions Si of A26[ satisfy ||5,||^i(i2) ^ '^"'^ 

Proof. Follows by a direct calculation. See Lemma lA3] for details. □ 



4 Stability of Galerkin discretizations 
4.1 Abstract results 

We consider the model problem (|9) and a sequence (yAf)ArGN C H^{Q) of finite- 
dimensional spaces. Furthermore, we assume that (Vw)a7gn is such that for every 
V £ H^{Q) we have limAr^ooinf,._^ev'A. Ih'^ ^'wllz/ifij) =0- The conforming approxi- 
mations Miv to the solution u of (|9]l are then defined by: 

Find MA, e Va, s.t. B{un,v)^1{v) Vv e Va,. (37) 

Since the sesquilinear form B satisfies a Garding inequality, general functional an- 
alytic argument show that asymptotically, the discrete problem (l37i has a unique 
solution UN and are quasi-optimal (see, e.g., I731 Thm. 4.2.9], 1741 ). More precisely, 
there exist A^o > and C > such that for all N>No 

\\u-UN\\i.k.n <C inf ||« - v|| 1.^,^2. (38) 

veVN 

This general functional analytic argument does not give any indication of how C 
and A^o depend on discretization parameters and the wavenumber k. Inspection of 
the arguments reveals that it is the approximation properties of the spaces Va, for the 
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approximation of the solution of certain adjoint problems that leads to the quasi- 
optimality result ( [38] l. For the reader's convenience, we repeat the argument, which 
has been used previously in, e.g., ||6][T3]|56]|6Tl|62l|72l|74l and is often attributed to 
Schatz, 1741 : 

Lemma 4.1 ( ll62l Thm. 3.2]). Let Q. dW' be a bounded Lipschitz domain and B 
be defined in ( 1771) . Denote by S* : L?{£2) — > //'(i2) the solution operator for the 
problem 

Findu* eH\n) s.t. B{v,u*) = {\'J)ii(a) &H\n). (39) 

Define the adjoint approximation property T7(Viv) by 

,(y.):= sup mfM!(^pfei. 

If, for the continuity constant Cb of M2^ . the space Vn satisfies 

2CB^77(y^)<l, (40) 
then the solution un oft[37\l exists and satisfies 

\\u-UN\\i,k.n <2Cb inf |1m - v|| (41) 

veVf^ 

Proof. We will not show existence of un but restrict our attention on the quasi- 
optimality result PTt ; we refer to [54^ Thm. 3.9] for the demonstration that i4T[ in 
fact implies existence and uniqueness of u^. Letting e = u — un he the error, we 
start with an estimate for ||e||i2(^): Using the definition of the operator S* and the 
Galerkin orthogonality satisfied by e, we have for arbitrary v eVn 

\\^\\l2{Q) = (^'^)L2(i2) ^ B{e,S*e) ^ B{e,S*e ~v) <CB\\e\\i,k,n\\S*e -vWi^k^Q. 

Infimizing over all v G Vn yields with the adjoint approximation property 77 (Vn) 

\\e\\LHn)<CBri{VN)\\e\\uk,n. 

The Garding inequality and the Galerkin orthogonality yield for arbitrary v e Vn: 

lkllUi2 = ReB(e,e) +2A:2||e||2,^^^ = ReB(e,M - v) + 2^2||^j|2^^^^ 

< 1 1 e 1 1 1 1 1 M - V 1 1 1 xn + {CBkri{VN)f \\e\\ ■ 

The assumption Cg^T] (Vn) < 1/2 allows us to rearrange this bound to get ||e|| i j^ q < 
2Cb\\u — v\\ i^k.n- Since v G Vn is arbitrary, we arrive at (1411 . □ 



Lemma \4m informs us that the convergence analysis for the Galerkin discretization 
of ^ can be reduced to the study of the adjoint approximation property J7(V)v), 
which is purely a question of approximation theory. In the context of piecewise 
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polynomial approximation spaces this requires a good regularity theory for the 
operator S*. The strong form of the equation satisfied by m* := S*f is 

-Au-k^u^f in 12, d„u*-iku'' = ondQ, (42) 

which is again a Helmholtz problem of the type considered in Section |3] More for- 
mally, with the solution operator S of Section [3] we have S*f ~ S{f,0), where an 
overbar denotes complex conjugation. Thus, the regularity theory of Section |3] is 
applicable. 



4.2 Stability of hp-F EM 

The estimates of Theorem [32] suggest that the effect of the corner singularities is 
essentially restricted to an (9 (1/fc) -neighborhood of the vertices. This motivates us 
to consider meshes that are refined in a small neighborhood of the vertices. To fix 
ideas, we restrict our attention to meshes i^/^" that are obtained in the following 
way: First, a quasi-uniform triangulation ,3^, with mesh size li is selected. Then, 
the elements abutting the vertices Ay, 7 = 1, ... ,7, are refined further with a mesh 
that is geometrically graded towards these vertices. These geometric meshes have L 
layers and use a grading factor a G (0, 1) (see 177] Sec. 4.4.1] for a precise formal 
definition). Furthermore, for any regular, shape-regular mesh we define 

SP{^):={ueH\n): uIk^^p VKe£^}, (43) 

where denotes the space of polynomials of degree p. We now show that on 
the geometric meshes .^f^", stability of the FEM is ensured if the mesh size li and 
the polynomial degree p satisfy the scale resolution condition (|6]l and, additionally, 
L = 0{p) layers of geometric refinement are used near the vertices: 

Theorem 4.2 (quasi-optimaUty of /z/j-FEM). Lef J'^f^' denote the geometric meshes 
on the polygon Q. <Z^ as described above. Fix C3 > 0. Then there are constants c\, 
C2 > depending solely on Q and the shape-regularity of the mesh such that 

the following is true: Ifh, p, and L satisfy the conditions 

kh 

— < ci and p>C2logk and L>cjp (44) 
P 

then the hp-FEM based on the space has a unique solution un G S''{3^i^"') 

and 

||"-Mw||i.i,r2 < 2Cb inf \\u-v\\x,k^a. (45) 

Proof. By Lemma|4T] we have to estimate ^77(Viv) with Vn = S'^.'T'if"'). Recalling 
the definition of r\ {Vn) we let f £L^{£2) and observe that we can decompose S*f = 
Ufj2 + Uj^, where m^2 and Uj^ satisfy the bounds 
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«//2|l//2(r2) < C||/!|^2(^), 



1)1+2 1 



V« e No. 



Piecewise polynomial approximation on =3^/^^'' as discussed in ||62| Prop. 5.6] gives 
under the assumptions kh/p < C and L > cj,p: (inspection of the proof of 162] 
Prop. 5.6] shows that only bounds on the derivatives of order > 2 are needed): 



inf ||M^2-V||ll,i2 < C-||/||i2(i2), 



inf ||m.^- v||i.<.,i2 < C 

veVN 



/ kh Y 



{Cso,ik)+l)\\fh2 



L^n) 



where j8„K„ = max,=i / j5j < 1, and C, c, b > are constants independent of h, p, 

and k. From this, we can easily infer 



kh 

kr]{VN)<C{j+k{Qo,{k) + l) 



{khf-^""''e''''''-''P + 



/ kh Y 



Noting that Theorem l2.4| gives Csoi{k) = 0{k^^^), and selecting ci sufficiently small 
as well as C2 sufficient large allows us to make A;77(Viv) so small that the condition 
(l40l i in Lemma|4T|is satisfied. □ 

Corollary 4.3 (exponential convergence on geometric meshes). Let f be analytic 
on Q and g be piecewise analytic, i.e., f, g satisfy ( I52I ). Given C3 > 0, there exist 
ci, C2 > such that under the scale resolution conditions Il44\l of Theorem \4.2\ the 
finite element approximation un G exists, and there are constants C, b >0 

independent ofk such that the error u — u^ satisfies 

||«-«A'||i,i-.i2 <Ce'''''. 



Proof. In view of Theorem 14.21 estimating \\u — unWi /^q is purely a question of 
approximability for ci sufficiently small and C2 sufficiently large. Lemma 1378] gives 
that the solution u = S{f,g) satisfies the bounds given there and, as in the proof of 
Theorem 14.21 we conclude from ||62] Prop. 5.6] (more precisely, this follows from 
its proof) 



inf ||m.c/- v||i,i,i2 < C 
veVN 



(kh^-P-'^'e^"-"' + ( —Y 



(C,„,(^) + l)(C/ + q). 



Theorem |2.4| asserts Csoi{k) = 0{k^l'^), which implies the result by suitably adjust- 
ing ci and C2 if necessary. □ 

Remark 4.4. 1. The problem size = dim5''(^ff ) is A? = 0{{L + h-^)p^). The 
particular choice of L ^ CT,p layers of geometric refinement, approximation 
order p — C2log^, and mesh size h = cip/k in Theorem 14.21 ensures quasi- 
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optimality of the hp-FEM with problem size i.e., quasi-optimality 

can be achieved with a fixed number of degrees of freedom per wavelength. 

2. The sparsity pattern of the system matrix is that of the classical hp-FEM, i.e., 
each row/column has 0{p^) non-zero entries. Noting that the scale resolution 
conditions (|6]l, (l44l i require p = 0{\ogk), we see that the number of non-zero 
entries entries per row/column is not independent of k. It is worth relating this 
observation to |12|. It is shown there for a model problem in 2D that no 9 point 
stencil can be found that leads to a pollution-free method. 

3. Any space Vn that contains where h, p, and L satisfy the scale reso- 
lution condition ( l44b also features quasi-optimality. 

4. The factor 2 on the right-hand side of (l45T l is arbitrary and can be replaced by 
any number greater than 1 . 

5. The stability analysis of Theorem |4.2| requires quite a significant mesh refine- 
ment near the vertices, namely, p. It is not clear whether this is an artifact 
of the proof. For a more careful numerical analysis of this issue, more detailed 
information about the stability properties of the solution operator S is needed, 
e.g., estimates for \\S{f,g)\\i_,,^Bu,(Aj)- 



4.3 Numerical examples: hp-FEM 

All calculations reported in this section are performed with the hp-FEM code NET- 
GEN/NGSOLVE by J. Schoberl, II75II76I . 

Example 4.5. In this 2D analog of Example 11.11 we consider the model problem 
(|9|l with exact solution being a plane wave e'('^i^+*2>'), where ki = —^2 = and 
k e {4,40, 100,400}. For fixed p e {1,2,3}, we show in Fig. |2] the performance 
of the /i-version FEM for p E {1,2,3} on quasi-uniform meshes by displaying the 
relative error in the //'-seminorm versus the number of degrees of freedom per 
wavelength. We observe that higher order methods are less prone to pollution. We 
note that the meshes are quasi-uniform, i.e., no geometric mesh refinement near the 
vertices is performed in contrast to the requirements of Theorem l4.2l ■ 

Example 4.6. On the L-shaped domain X2 = (-l,l)2\(0,l)x(-l,0) with F being 
the union of the two edges meeting at (0,0), we consider 

— Au — k^u^O ini2, d„u = onF, d„u — iku = g on dQ\r, (46) 

where the Robin data g are such that the exact solution is u{x,y) = e'(*^i-''+'^2>') 
with ki = —ki = -^^1 ™d ^ e {10, 100, 1000}. We consider two kinds of meshes, 
namely, quasi-uniform meshes ^1, with mesh size h such that kh w 4 and meshes 
(^geo fjj^j- ^j.g geometrically refined near the origin. The meshes ^s^o ^j.g derived 
from the quasi-uniform mesh 5/, by introducing a geometric grading on the ele- 
ments abutting the origin; the grading factor is a = 0.125 and the number of refine- 
ment levels is L = 10. Fig. [3] shows the relative errors in the //^-seminorm for the 
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p=1, plane wave solution, unit square p=2, plane wave solution, unit square 




DOF per wavelength DOF per wavelength 




DOF per wavelength polynomial degree p 



Fig. 2 Top: /i-FEM with p = 1 (left) and p = 2 (right) as described in Example 14.5! Bottom left: 
/i-FEM with /) = 3 as described in Example |4.5| Bottom right: p-FEM for singular solution on 
quasi-uniform mesh as described in Example 14. 7 1 



p-version of the FEM where for fixed mesh the approximation order p ranges from 
1 to 10. It is particularly noteworthy that the refinement near the origin has hardly 
any effect on the convergence behavior of the FEM; this is quite in contrast to the 
stability result Theorem |4.2| which requires geometric refinement near all vertices 
ofQ. m 

Example 4.7. The geometry and the boundary conditions are as described in Exam- 
ple |4j6| The data g are selected such that the exact solution is u = J2/3{kr)cos 
where (r, (j?) denote polar coordinates and is a first kind Bessel function, k G 
{1,10,20, 100,200}. Our calculations are p-FEMs with p G {1, . . . , 10} on the qua- 
siuniform mesh 5^, described in Example 14.7 1 The results are displayed in the bot- 
tom right part of Fig. |2] The numerics illustrate that the singularity at the origin 
is rather weak: we observe that the asymptotic algebraic convergence behavior is 
\u — MAf|/fi(^) ~ Ckp^'^^^\u\un^Q-^, where the constant Q depends favorably on k. ■ 
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plane wave solution, L = 0, a = 0.125 plane wave solution, L = 10, o = 0.125 
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Fig. 3 p-FEM for plane wave solution as described in Example |4.6l Left: quasiuniform mesh 
with kh Ri 4. Right: Mesh obtained from .% by strong geometric refinement near origin. 



4.4 Stability of Partition of Unity Method/Generalized FEM 



The abstract stability result of Lemma l^Tl onlv assumes certain approximation prop- 
erties of the spaces Vn- Particularly in an "/i-version" setting, even non-polynomial, 
operator-adapted spaces may have sufficient approximation properties to ensure the 
important condition (|40] | for stability. We illustrate this effect for the PUM/gFEM, 
Il56ll60l with local approximation spaces consisting of systems of plane waves or 
generalized harmonic polynomials (see Section |5] below) and the classical FEM 
shape functions as the partition of unity. The key observation is that for li sufficiently 
small, the resulting space has approximation properties similar to the classical (low 
order) FEM space: 

Lemma 4.8. Let be a shape-regular triangulation of the polygon £2 C M?. Let 
h be its mesh size; let ^''^ nodes of the triangulation and be the 

piecewise linear hat functions associated with the nodes (x, )^j. Assume kh < C\ 
for some Ci > 0. Let V"'"""''' be either the space V^f^p with p>0 (see below) 
or the space Wp^ with p > 2 (see ([50i below). Define, for each i ~ 1, . . . ,M, the 
local approximation V,- by V,- :— span{b{x — Xi) : b e ymastery j^i^gji f/jg space Vn :— 
^h=i ^i^i ^'^^ the following approximation property: There exists C > depending 
only on the shape regularity of , the constant C\, and v'"'"'""' such that 



inf \\u-v\\LHn 
veVN ^ 



-h\\u — v\ 



<C 



Proof. We first show that each local space V, has an element Xj/j E V, with 



/2||Vv/,| 



\^-Wi\ 



L°°((a,) 



^ Capp{kh 



(47) 



for some Capp > independent of ; and h; here, d), = supp^,- has diameter 0{h). 
It suffices to show ( |47] | for the set y"""^"'''. For the space of generalized har- 
monic polynomials, this follows from Jo{kr) = 1 + 0{{kr)~), and for the space 
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of plane waves, Taylor expansion shows that for p ~ 2m (m e N) the function 
2 [e'*^* + e^'^'^] = 1 +0{{kx)^) has the desired property whereas for odd p = 2m + 1 
(m € N) the observation 



2c 



,!<:(»,„ A I i<r(a„,+ i -.r 



1\ 27Zm 



can be utilized to construct i/,-. We recall (see, e.g., ifTSl Thm. 4.8.7]) that for each 
u e H^{£2) there is a function w — ^^(^{^^(Pi G 5' (3^,) with 

M\L^(n)<C\\u\\[2^Q), Ikll//i(i2) (48a) 

\W~w\\LHn)<Ch^h\\H^{n)^ \W-w\\H>{n) <Ch\\u\\H2^Qy (48b) 

Upon setting v := LiM^K^i) Vif! G ^n, we get in view of = 1 for the error 
u~v=u — 'E,iWi(xi)WiVi = (m — w) + Li' (Piw{xi){l — I/;). The estimates ( l48T l imply 
Ik' ^ '^llL2(i2) — '^Il//i(i2) ^ C/i^llMll^a^ij). For the sum, we have 

\\Y,Wi{xi){l-\l/,)(p,\\^2(n) < C(fc/l)^||w||^2(i2) <C(^/j)^||M|lz,2(i2), 
/ 

/i|lV^w,(x,-)(l - r,)<P-llL2(r2) < CikhfM^^nn) < Cikhf\\u\\i^2^a), 

i 

which concludes the proof. □ 

Remark 4.9. The approximation result of Lemma 14.81 can generalized in various 
directions. For example, interpolation arguments allow one to construct, for v e 
/fi+^(i2) with 9 e (0, 1) an approximation Vap,, G Vn such that \\v — Vapp\\L2(n) + 
h\\v-Vapp\\HHn]<Ci (/i'+'' + (^/i)2/i9)|M|;,i+e(^2)+C2(fc/i)2||«||^2(^). (We refer 
the reader to the Appendix for the proof of this results.) Furthermore, the approxi- 
mation result of Lemma |4~8] can be localized, which is of interest if 5^ is not quasi- 
uniform. ■ 



Lemma l478l shows that the space Vn, which is derived from solutions of the homoge- 
neous Helmholtz equation, nevertheless has some approximation power for arbitrary 
functions with some Sobolev regularity. Hence, the condition (|40] | can be met for 
sufficiently small mesh sizes: 



Corollary 4.10 ( 11561 Prop. 8.2.7J). Assume the hypotheses of Lemma \4.8\ in par- 
ticular, let the space be constructed from systems of plane waves or generalized 
harmonic polynomials. Assume additionally that Q is a convex polygon. Then there 
exists C > independent ofk such that for k^h < C the Galerkin method for ([9l) with 
f = is quasi-optimal, i.e., the solution un G Vn of 071) exists and satisfies 

\\u-UN\\i.k.n <2Cb inf ||m - v|| n-.^^. 

veVN 



Proof. In view of Lemma BTTl we have to estimate Tj (VAf). To that end, we consider 
(|9]l with / e L^{£2) andg = 0. In view of the convexity of i2, we have Csoi{k) = (9(1) 
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and elliptic regularity then yields for the solution u of (|9]l 

\W\\i.k.n+k''^\\u\\H2{n) <C\\f\\[2(^ay 
This allows us to conclude with Lemma l4~8] that 

{kh^+ll)\\u\\f,2^^) + {k{kh)^+k^h)\\u\\^2^^^ 
< C((^/!)2+A:/0)||/||^2(^) <C^/l(l+W7)||/||^,(^). 

Hence, A;T7(V/v) can be made sufficiently small if k^h is sufficiently small. We point 
out that convexity of Q is assumed for convenience — under more stringent condi- 
tions on the mesh size li, quasioptimality holds for general polygons. □ 



inf I 

VGVm 



< C 



5 Approximation with plane, cylindrical, and spherical waves 

Systems of functions that are solutions of a (homogeneous) differential equation are 
often called "Trefftz systems". Prominent examples in the context of the Helmholtz 
equation are, in the two-dimensional setting, "generalized harmonic polynomials" 
and systems of plane waves given by 

V^HP ■■= span{y„(^r)e'"'P : -p<n< «}, (49) 
W/n, : = spanje'*-""" : « = 0, . . . , p - 1 } , (0„ = (cos ^ , sin ^ ) ; (50) 



here, J„ is a first kind Bessel function, the functions in VA^^ are described in polar 



coordinates and the functions of Wp^, in Cartesian coordinates. We point out that 



'GHP 

analogous systems can be developed in 3D. These functions are solutions of the ho- 
mogeneous Helmholtz equation. For the approximation of a function u that satisfies 
the homogeneous Helmholtz equation on a domain D. C M^, one may study the "/?- 
version", i.e., study how well u can be approximated from the spaces V^^p or 
as /7 ^ oo; alternatively, one may study the "/i-version", in which, for fixed p, the 
approximation propeities of the spaces Vq^p or W^jy are expressed in terms of the 
diameter h = diami2 of a domain under consideration. In the way of illustration, we 
present two types of results: 

Lemma 5.1 ( 11561 ). Let £2 C M? be a simply connected domain and £2' CC £2 be 

a compact subset. Let u solve —Au — k^u = on Q. Then there exist constants C, 
b > (possibly depending on k) such that for all p > 2: 



inf ||«-v||^,(^,)<Ce-^ inf |1« - v||^,(^,) < Ce-^'^mp. 



Proof. See, e.g., ||56l or EH Thm. 5.3]. □ 
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Remark 5.2. Analogs of Lemma lSTI hold if u has only some finite Sobolev regular- 
ity. Then, the convergence rates are algebraic, ll56l . Il58] Thm. 5.4], 1421 . ■ 

The approximation properties of the spaces V^^p and Wp^, can be also be studied in 
an /j-version setting: 

Proposition 5.3 ( ll42l Thm. 3.2.2]). Let Q cM? be a domain with diameter h and 
inscribed circle of radius ph. Let p = 2jX + \. Assume kh < C\. Then there exist 
Cp> (depending only onC\, p > 0, m, and p) and v G Wp^^^ such that 

h- v\\j,k,n,E < Cph^''^^\\u\\^,+l.k.nx, 0<j<fi + l, 

Several comments concerning Proposition |53]are in order: 

1. The constant Cp in Proposition l5.3l depends favorably on p, and its dependence 
on p can be found in 142] Thm. 3.2.3]. 

2. Proposition 15.31 is formulated for the space Wp^, of plane waves — analogous 
results are vahd for generalized harmonic polynomials, see ll42l Thm. 2.2.1] for 
both the h and /jp-version. 

3. Proposition [53] is formulated for the two-dimensional case. Similar results are 
available in 3D, l42l . 

4. The approximation properties of plane waves in terms of the element size have 
previously been studied in slightly different norms in ||20l . 

Remark 5.4. Plane waves and generalized harmonic polynomials represent by no 
means the only operator adapted system used in practice. Especially for polygo- 
nal domains, the functions J„a{kr) sm{an(p), n G N, or J„a{kr) sm{an(p), n G No, 
for suitable a can combine good approximation properties with the option to re- 
ahze homogeneous boundary conditions, lT4l . Further possibilities include linear 
combinations of fundamental solutions or, more generally, discretizations of layer 
potentials. We refer to |[T4l for a concrete example. ■ 



6 Stability of Least Squares and DG methods 

Discrete stability in Section |4] is obtained from stability of the continuous prob- 
lem by a perturbation argument. This approach does not seem to work very well 
if one aims at using approximation spaces that have special features linked to the 
differential equation under consideration. The reason can be seen from the proof of 
Lemma \4~T\ The adjoint approximation property T7(Va?) (which needs to be small) 
measures how well certain solutions to the /«homogeneous equation can be approx- 
imated from the test space. If the ansatz space is based on solutions of the homoge- 
neous equation, then its capabilities to approximate solutions of the inhomogeneous 
equation are clearly limited. In an /i-version, the situation is not as severe as we have 
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just seen in Section |4~4l for the PUM/gFEM. In a pure /9- version setting, however, 
the techniques of Section |4!4l do not seem applicable. 

An option is to leave the setting of Galerkin methods and to work with formu- 
lations with built-in stability properties. Such approaches can often be understood 
as minimizing some residual norm, which then provides automatically stability and 
quasi-optimality (in this residual norm). We will illustrate this procedure here by 
two examples, namely. Least Squares methods and DG-methods. Our presentation 
will highlight an issue stemming from this approach, namely, the fact that error es- 
timates in this residual norm do not easily lead to error estimates in more classical 
norms such as the L^(i2)-norm. 



6.1 Some notation for spaces of piecewise smooth functions 

Let 5^ be a regular, shape-regular triangulation of the polygon 12 C K^. We decom- 
pose the set of edges as ^ = S'jOSb, where Si is the set of edges in Q. and S'b 
consists of the edges on dO, . For functions m : 12 — > R and (7 : 12 — > that are 
smooth on the elements iiT € =5^, we define the jumps and averages as it is customary 
in DG-settings: 

• For e € let and be the two elements sharing e and denote by n+ and 

the normal vectors on e pointing out of and . Correspondingly, we let 
M+, and (7+ and CJ^ be traces on e of u and a from and . We set: 

iu}}\e \ («+ + U-) , ia^V \ + c^") , 
\u^\e := M+n+ + M"n", [[(y]]\e := CT+ • n+ + f7" ■ n". 

• For boundary edges e £ S'b we define 

IcT^-le := [[m]]|<. := mI^ii 

With this notation, one can conveniently rearrange certain sums over edges; 

Lemma 6.1 ("DG magic formula")- Let v ; i2 — > K and a : £2 be piecewise 

smooth on the triangulation Then: 

£ / vcT • n = / M . {{dS + / M ■ M + i M ■ i<y^. 

where jg^ and Jg^ are shorthand notations for the sums of integrals over all edges 
in Si and Sb- 

Finally, for piecewise smooth functions, V/, denotes the piecewise defined gradient. 
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6.2 Stability of least squares methods 

Although Least Squares methods could be based on any space of approximation 
spaces, we will concentrate here on the approximation by piecewise solutions of 
the homogeneous Helmholtz equation. With varying focus, this is the setting of 
II 1 4115 3 Il65 Il70ll80ll and references therein. We illustrate the procedure for the model 
problem (|9|l with / = 0. The approximation space has the form 

Vn^{u(^L^{Q):u\k(^Vn,k V/Te^}, (51) 

where the spaces Vn^k are spaces of solutions of the homogeneous Helmholtz equa- 
tion, e.g., systems of plane waves. For each edge e G (a", we select weights wi g, 
W2.e > and define the functional J : Vn ^ M.hy 

7(v):= £ wlMy]\\l2^,)+^ie\\[d,M\\l2^,)+ E ^ie\\8-idnv + ikv)\\l2^^y 

eeSj eeSB 

here \y\\e '■= [[v]] \e and [d„v] \e'~ [[V/,v]] \e represent the jumps of v and d„\' across the 
edge e. If the exact solution m of (|9]l is sufficiently regular, then it is a minimizer of / 
with J{u) = 0. In a Least Squares method, J is minimizer over a finite dimensional 
space Vn of the form ( |5TT i. Its variational form reads: 

find MA, e Va, s.t.(MA,,v)/,A< = 21 (g,^„v + ifcy)^2(g) Vv e V^v, (52) 

where 

{u,v)j,N :== 

The positive semidefinite sesquilinear form {■,-)j,n induces in fact a norm on Vn: 
To see the definiteness of (•, ■)j.n, we note that v E Vn and J{v) = implies that v 
is in C'(f2) and elementwise a solution of the homogeneous Helmholtz equation. 
Thus, it is a classical solution of the Helmholtz equation on Q and satisfies d„v + 
ikv = on d£2 . The uniqueness assertion for (|9]l with / = and ^ = worked out 
in Example 12.11 then implies v = 0. Therefore, the minimization problem (l52t is 
well-defined. If the solution u of ^ satisfies ii G H^/^+^ (i2 ) for some £ > 0, then 
J{u) — 0, and we get quasi-optimality of the Least Squares method in the norm 

II -lU = 

IIM-Mwll/Af ^J{u-un) =J{un) = min7(v) = ||m- vjl^Af. (53) 

vev^ 

We mention here that estimates for this minimum can be obtained from (local) es- 
timates in classical Sobolev norms as given in Section |5] using appropriate trace 
estimates. Turning estimates for \\u — un\\j,n — -/(mw)'''^ into estimates in terms of 
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more familiar norms such as ||m — MwH^aj^j is not straight forward. It may be ex- 
pected that the norm of the solution operator of the continuous problem appears 
again; the next result, which is closely related to |fT9]l42l|43ll63l . illustrates the kind 
of result one can obtain, in particular in a p-version setting: 

Lemma 6.2 ( [|65t Thm. 3.1]). Let £2 a be a polygon. Let w\^e~k and W2,e = 1 
for all edges and g £ L^{dQ). Let un G Vn be the minimizer of J, where V)v, given 
by ( 157 H . consists of elementwise solutions of the homogeneous Helmholtz equation. 



(i) If Q is convex, then ||m — MAr||^2(^-) < CA: 

(ii) If Q is not convex, then 



-1 



{{kh)-^ + ikhy)j{uN). 



||M-MA'|li2(^2- 



< 



Ck 



-1 



(kh)-^ + |l +min{l,^/j}-2fi™.}] (Qoiik) + 1)^J{un), 



where Csoi{k) is defined in ( |25| ) and satisfies Csoi{k) = 0{k^^^) by Theorem \2.4\ 
The parameter ^,nax ^ can be selected arbitrarily to satisfy the condition 
Pmax > 1 ~ min,- where the (Oi are the interior angles of the polygon. 



Proof. The result Q is essentially the statement of ll65l Thm. 3.1] in a refined form 
as given in 1431 Lemma 3.7]. The statement is a slightly modification of and 
we restrict our presentation to that case. The key idea is to obtain L'^(^2)-bounds by 
a duality argument and use the fact that u — un solves the homogeneous Helmholtz 
equation elementwise. More precisely, given ^ G L^(i2) we let v G //' (i2) solve the 
adjoint problem 

—Av — k^v = (p in 12, d„v — ikv = ondQ. 



By Corollary 13.31 the function y is in a weighted // -space with 

I|v||i,^,i2 +^-'ll^o,^,v2v||^2(^) < C(C,„,W + l)\\(p\\L2iny (54) 

Inspection of the arguments underlying the proof of Corollarv l3.3| shows that the ex- 
ponents Pj G [0, 1) stem from the regularity theory for the Laplacian with Neumann 
boundary conditions. Hence, in fact j3/ e [0, 1/2) so that V/,v has an -trace on all 
edges of the triangulation (cf. Lemma IaT^I i. For each /iT e 5^ we then have 



w\ 



< C 



IVwII?, 



L-{dK) 



< C 



'llHlt2(^)+/l|w|?„ 



h-'\w\l, 
';,-l|..,|2 



HHK) 



VweH^K), 
\fweH^{K), 



Vw 



eH^-^K), (57) 



(55) 
(56) 



where, in the last estimate we assume that the origin is at one corner of K and r 
denotes the distance from that corner These estimates are obtained with the aid of 
scaling arguments, the multiplicative trace inequality (see ifTSl Prop. 1.6.3]), and. 
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in the case of ( |57] ) additionally Lemma pV. 21 From ab = max{a,b}mm{a,b} (for 
a,b> 0), we get rP = k-P{rkf = k'^ mm{\, rk}^ max{l,ryt}^. Hence, if ^T^omer 
denotes the set of elements that abut on the corners of i2, we get 



Ke,% 



corner 



where B/, = ^K(ifircor„„^- Noting that the elements K £ ,!^ \ ,%orner are at least 0{h) 
away from the corners allows us to estimate with ( 



Hence, we have the two bounds 



0,/3,A: 'llL2(r2\B,,) 



L l|Vv|j2, < C/r'|v|^, +C/.'min{l,^M-2A™-'||0^^^v2v||2, 

L Ml2iSK)<Ch-'\\v\\l2^^^+h\v\l,^^y 

Therefore, recalling w\,e = k and W2.e = 1, we obtain from these estimates and the 
a priori estimate ( l54b the bound 



eeS" 



< Ck- 



— +kli{l+ min{ 1 , khy^f^'""' ] 
kh L J 



{CUk) + ^fMn,nv (58) 



The estimate (ISSl l can be used to bound \ {u — UAf, ^)^2(^) |: Writing the integral as 
a sum over elements, integrating by parts twice and using that u — un solves the 
homogeneous Helmholtz equation elementwise yields 

{u-UN,(p)L2(n) = L {u-UN,-Av-k^v)^HK) 
The "DG magic formulas" of Lemma|6T|produce 

eG<5/ e^Sg 



For interior edges e € if/ we have [[v]] = and [[V\']] = as well as [[m]] = and 
[[Vm]] = 0; on boundary edges e G (og we have with the boundary conditions satisfied 
by M and v (i.e., (3„m + iku = g and (9„ v — ikv ~ 0) 



On stability of discretizations of the Helmholtz equation (extended version) 



35 



These observations lead to 



-2ii fr..Ti 112 , ...-2ii 

V 

where, in the last step, we employed the Cauchy-Schwarz inequality for sums. From 
we therefore get 

(M-MjV,(p)^2(^)| ^ 



\m\L2{n) 

Cy/j{uN)k-^'^ Vkh)-"^!^ + {khYl^mm{l,kh}-^""A {Qoi{k) + 1). 
Since 9 g (12 ) is arbitrary, we get the result. □ 

Remark 6.3. Lemma l6^ assumes quasi-uniform meshes and the weights w\ e, W2,e 
do not take the edge length into account. This limits somewhat it applicability in 
an /z-version context. However, the result is very suitable for a /j-version setting. 
We point out that in a case where the p-version features only algebraic rates of 
convergence, one would have to give the parameters wi^, W2,e a p-dependentrelative 
weight as opposed to the situation studied in Lemma l672l ■ 



6.3 Stability of plane wave DG and UWVF 

The framework of Discontinuous Galerkin (DG) methods permits another way of 
deriving numerical schemes that are inherently stable. In a classical, piecewise poly- 
nomial setting, this is pursued in ||33H35l ; related work is in 1641 . Here, we concen- 
trate on a setting where the ansatz functions satisfy the homogeneous Helmholtz 
equation. In particular, we study the plane wave DG method, [|36l'43"631, and the 
closely related Ultra Weak Variational Formulation (UWVF), 119-21 , 46, 55J. We 
point out that the UWVF can be derived in different way. Here, we follow II19II36I 
in viewing it as a special DG method. 

Our model problem (|9) can be reformulated as a first order system by introducing 
the flux a: = {l/ik)Vu: 



ika = Vu in £2, 
i^M — V • (7 = in i2 , 
ika-n + iku = g ondQ. 



(59a) 
(59b) 
(59c) 
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For a mesh ST , the weak elementwise formulation of ( |59a| i. ( |59bl l is for every K e 
f \ka-T+ [ mV-t- / MT-n = Vt e //(div,ii:), 

K Jk JdK 

ikuv+ [ a-Vv- [ c7-nv = \fveH\K), 
K Jk JdK 

where H{di\,K) = {« e L^^^) . ^^^^^^ ^ and n is the outward pointing nor- 

mal vector. Replacing the spaces H^{K) and H{di\,K) by finite-dimensional sub- 
sets Vn,k C H^{K) and Ln.k C //(div,/r) and, additionally, imposing a coupling 
between neighboring elements by replacing the multivalued traces u and a on the 
element edges by single- valued numerical fluxes un, On to be specified below, leads 
to the problem: Find {un, On) G Vn,k x ^n.k such that 

ikGN-^+ Un'^-X—1 MArT-n = Vt E En K, 

k Jk JdK 



L 



ikuNV+ I (7Af-Vv— / c7Af-nv = '^veVnk- 
K Jk JdK 

The variable On can be eliminated by making the assumption that Wn.k C T,n.k 
for all A" G ^ and then selecting the test function T = Vv on each element. This 
yields after an integration by parts: 

/ yuNyv — k^UNV— / {iiN -UN)d„v -ikON -nv = Q V/iT e 5^. (60) 
Jk JdK 

Since Vn = {u E L^{£2): u\k G Vn.k'^K e consists of discontinuous functions 
without any interelement continuity imposed across the element edges, ( l60l l is 
equivalent to the sum over the elements: Find un EVn such that for all v E Vn 

V / Vun ■'Vv — k^UNV+ / (i?Ar — ma') • n — / ikON -i^y = 0. (61) 

Jk JdK JdK 

This formulation is now completed by specifying the fluxes un and On, which at 
the same time takes care of the boundary condition (|59c) : 

• For interior edges e £ S"/ 

On = ;7{{V/,Miv}} - a [[un]] , UN = imj -Prr [[V/i«A']] • (62a) 
ik ik 

• For boundary edges e E S'g 

Y ^ 

On rr'^hUN rr- C^hm + ikuNii - gn) . (62b) 

ik ik 

UN = UN (V/,M • n + ikuN — g) ■ (62c) 

ik 
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Different choices of the parameters a, j5, 5 lead to different methods analyzed in 
the Uterature. For example: 

1 . a = j3 = 5 = 1/2: this is the UWVF as analyzed in |fT9ti2l1l46ll55l if the spaces 
Vn,k consist of a space Wf^, of plane waves. 

2. a = 0{p/{kh\ogp)), P =0{{khlogp)/p), 5 = 0((/fe/zlogp)//7): this choice 
is introduced and advocated in 04311631 in conjunction with Vn,k = ^/vr- 

With these choices of fluxes, the formulation i6l[ takes the form 

Find M^f e V^f s.t. Afj{ufj,v) = l{v) Vv € V/v, (63) 
where the sesquiUnear form and the linear form I are given by 

An{u, v) = j ■ Vhv-k^iiv- jMV^hV^ - f i^huW^ - 1 Sud„v- [ 5dnuv 

J £2 J S^i J Si J Sg J Sg 

-It I i3[[V/,M]][[V/,y]]-i / 5d„ud„v + ik [ «[[«]][[?]] (l-5)«v (64) 



So far, the choice of the spaces Vn.k is arbitrary. If the approximation spaces Vn.k 
(more precisely: the test spaces) consist of piecewise solutions of the homogeneous 
Helmholtz equation, then a further integration by parts is possible to eliminate all 
volume contributions uiAn. Indeed, Lemma l6?T] produces 



so that Aa? simplifies to 



dK Jg, Js, 



An{u,v)= I M[[V,v]]+ii/ j3[[V,«]][[V,,v]]_ / |v,4[[v]]+i^/ aH[[v]] 

J Si 1^ JS/ J S'l J Sj 

(1 — 5)«(?„v + i- / dd„ud„v— [ 5d„uv + ik ( (1— 5)mv. 



k J St 

Next, we make the important observation that ImAAf induces a norm on the space 
ViV if a, j8 > and 5 e (0, 1). Indeed: 

1. a, j3 > and 5 G (0, 1) implies ImAAf(v, v) > Vv £ Vn bv inspection of (l64b. 

2. ImAAf(v, v) = and the fact that Vn consists of elementwise solutions of the ho- 
mogeneous Helmholtz equation implies as in the case of (•, •)yjv in Section [672 
that V e C'(i2) solves the homogeneous Helmholtz equation and d^v v = 
on dQ ; the uniqueness assertion of Example 12 . II then proves v = 0. 

This is at the basis of the convergence analysis. Introducing 
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\\u\\lc ■■= v/lmA^(«,«) = ^lli3'/'[[V/,«]]|l22(^^) + ||aV2[[«]]||2,^^^^ 
we can formulate coercivity and continuity results; 

Proposition 6.4 ( II19II43I ). Let Vn consist of piecewise solutions of the homoge- 
neous Helmholtz equation. Then \\ ■ \\j)q is a norm on Vn and for some C > de- 
pending solely on the choice of a, ji > 0, and 5 G (0, 1); 

lmAN{u,u) = ||m||dg ^« £ ^N, 

\An(u,v)\ < C||M||DG,+ l|t'l|DG Vm,v e Vn 

Let the solution of u of ^ (with f ~Q) satisfy u S H^^'^^^{Q) for some e > 0. 
Then, by consistency of An, the solution u^ G Viv of ( 1631 ) satisfies the following 
quasioptimality estimate for some C > independent of k: 

\\u-un\\dg <C inf ||«-vj|DG,+- (65) 

Several comments are in order: 

1 . The UWVF of 1201 featured quasi-optimality in a residual type norm. We recall 
that the UWVF is a DG method for the particular choice a = j3 = 5 = 1/2. 

2. When Vn consists (elementwise) of systems of plane waves or generalized har- 
monic polynomials, then the infimum in ( |65] | can be estimated using approxi- 
mation results on the elements by taking appropriate traces. This is worked out 
in detail in i42l l43] , l63l and earlier in an /i-version setting in 120} (see also [|T9ll ). 

3. The II • lloG-norm controls the error on the skeleton S only. The proof of 
Lemma 16.21 shows how error estimates in such norms can be used to obtain 
estimates for ||m — un\\iiiq\; we refer again to |fT9l where this worked out for 
the UWVF and to Il42ll43ll63l where the case of the plane wave DG is studied. 
As pointed out in Remark l631 quasi-uniformity of the underlying mesh ,^ is an 
important ingredient for the arguments of Lemma |6^ 

It is noteworthy that Proposition 16.41 does not make any assumptions on the mesh 
size h and the space Vn except that it consist of piecewise solutions of the homo- 
geneous Helmholtz equation. Optimal error estimates are possible in an /i-version 
setting, where the number of plane waves per element is kept fixed: 

Proposition 6.5 ( |l36l). Let Q be convex. Assume that Vn.k = Wpl^^^ (IX > I fixed) 
for all K ^ S/' . Assume that (X is of the form a ~ a/{kh) and that j3 > 0, 5 G (0, 1 /2). 
Then there exist ao, Co, C > (all independent of h and k) such that if &> ao and 
k^h < Co, then following error bound is true: 

\\u-un\\i,dg < C inf ||m - v|| i.og,+; 
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/ze'-e, 1 1 ■ 1 1 LDG an^ 1 1 J 1 ,OG. + fl'-e ven 1 1 V 1 1 LDG ■ = I/f e .3^ I ^ I Hi (^) + 1 1 ^' 1 1 l2 (/f ) + 
MIg MloG.^ := I^e^ MIi(k)+'^'Ml2{k) + Mdg..+- 

Proof. The proof follows by inspection of the procedure in ll36] Sec. 5] and is stated 
in 1631 Props. 4.2, 4.3]. The essential ingredients of the proof are: (a) inverse esti- 
mates for systems of plane waves that have been made in available in |[36l so that 
techniques of standard DG methods can be used to treat A^; (b) use of duality ar- 
guments as in Lemma 143] to treat the L^-norm of the error; (c) the fact that in an 
/i-version setting, plane waves have some approximation power for arbitrary func- 
tions in (this is analogous to Lemma |48] |. □ 



7 Remarks on ID 

The ID situation is rather special in that pollution can be completely eliminated; 
the underlying reason is that the space of solutions of the homogeneous Helmholtz 
equation is finite-dimensional (two dimensional, in fact). We illustrate this for the 
following model problem: 

-u"-k^u=f ini2 = (0,l), m(0)=0, u'{\)-iku{\)^Q. (66) 

Let 5^ be a mesh on Q. with nodes 0=:JCo<-*^i < ■ ■ ■ < = I. We assume that the 
mesh size h := max,(x,+i — x,) satisfies kli < n. For each node x,-, let y/i e H^{£2) 
be defined by the conditions 

Wiixj) = 5,j, i-Wi - k^Vi)\K = yxe-^r 

and let V^'" = span{ y/, : ;=!,... ,A^}. Thus, the functions Xj/j are piecewise solutions 
of the homogeneous Helmholtz equation. The Galerkin method based on V^'^' is: 

Find MA, e y^'" s.t. [ u'j^v' -k^UNV-\kuN{\)v{\)= ( fv VvSV^'". (67) 
Jn Jn 

The Galerkin method based on V'^'^' is nodally exact: 

Lemma 7.1. There exist constants C\, C2 > independent ofk such that the follow- 
ing is true for kh < Ci: 

(i) The functions are well-defined. 

( ii) The method A67[ is nodally exact. 

(Hi) For f eL?{£1) there holds Wu — UN\\\.k,n 1^C2{hk)\\f\\i2(^Qy 

Proof. Elementary considerations show that for kh < n, the functions xj/j are well- 
defined. 

The most interesting feature of Lemma ItTI is the nodal exactness. To that end, 
we note that the Green's function for (|66] | is 
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, , 1 fsin^jce'*^' 0<x<y 
G(x,y) = - < ., 
^ ^ k [sin/tye''^^ 3;<x< 1. 

Let e G //' (i2) satisfy e(0) =0 and the Galerkin orthogonality condition 



C(e,y):= / ev-k^ev-ike{l)v{l)=0 Vv G 
Jn 



opt 
N ■ 



(68) 



The key observation is that for each x;, / = 1, . . . ,A^, the function v,- G(-,x,) G V^'" 
since is a solution of the homogeneous Helmholtz equation on (0,jc,) U (x,, 1), it 
satisfies G(0,jc,) = 0. Furthermore, we have v'^{xn) — ik\'i{xi\i) = 0. Hence, we get 
from the Galerkin orthogonality (l6Ft by an integration by parts: 



= / e'v'i-k^eVi + ike{\)vi{\) 
Jn 

= / ev'i-k^evi+ / e'v'i-k^evi-ikeiiyviil) 

Jo Jxi 
pi 

= I e{-\'l -l?Vi)+e{xi)\'\{xi)-e{xi)v\{xi) 
Jo 

+ [\{-v'!-kh,)+e{l)x^,{l)-ike{l)v,{l) 

Jxj 

= e{xi)[v']{xi)+e{imi)-ike{l)vi{l) 
= e{x,M{x,y, 

here, we have employed the standard notation for the jump of a piecewise smooth 
function w: [w'] (x,) limj^x,- — limv^v.+ w'(jc). Since [v'] (x,) 7^ 0, we con- 
clude 

e(x,) = V/G {!,..., A^}. 

Hence, the FEM (|67] ) is nodally exact. 

The above argument also shows that any e G V^'" satisfying ( |68] | must satisfy 
e{xi) ~ for all ; G {1, • • • ,A^}- Hence, by the definition of V^'" as the span of the 
functions we conclude e = 0. Thus, the kernel of the linear systems described by 
([67]) is trivial. By the usual dimension argument, we have unique solvability. 

We have the a priori bound 

Mi.k.n<C\\f\y^n^ (69) 

for the solution u of ( |66] | (as for the model problem (|9]l, this can be shown using 
the test function v ~ xu'; an alternative proof based on the Green's function and the 
representation 

u{x) = / G{x,y)f{y)dy 
Jn 

is given in BTl Thm. 4.4]). If one denotes by (pj the classical piecewise linear hat 
function associated with node Xi, then one has by Taylor expansion 
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\\(P, - w.h-iK) < c{khK)\ \\i(Pi - wd'\\L-(K) < ck^zK, yK G sr, (70) 

where — diam/ir. The approximation properties follow easily from the nodal 
exactness. Specifically, denoting by lu the classical piecewise linear interpolant of 
u and by lu € V^'^' the nodal interpolant determined by lu(xi) ~ u(xi) for all / e 
{0, . . . ,A^}, we have the well-known estimate 



II — III 



Wk-a < C{kh' + h)\\ii"hHn) < Ckh{\ +kh)\\f\\^,^n^, 



where we used the differential equation and the bound (|69) to estimate ||«"||i2(^) < 

ll/llL2(r2) +^^lk'llL2(i2) ^ C^\\f\\L-(n)- Next, we estimate the difference /m — /m. The 
multiplicative trace inequality takes the form 

/%||w|l^(^) <c[||w||2,(^,+/,^||w||^2(^)||w'||i2(^)] yw€H\K). (71) 

Hence, the estimates dTOl i. ( ItTI i imply 



\I>^~M\ik.n<C k^{khKf{\ + {khKf)hK\\u\\ 



2 

Ke.sr 



< c{khY{\ + [khY) \k'\\u\\l.^^^+k'h%u'\\l^^^^ 

<C{khf{\ + {khff\\u\\l,^a- 
An appeal to (|69] | concludes the argument. □ 

Several comments are in order concerning the stability of the method: 

1. In the ID situation, the good stability properties of high order Galerkin FEM 
can alternatively be understood in light of Lemma ItTI Applying the Galerkin 
method to a classical high order method and then condensing out the degrees 
of freedom corresponding to internal shape functions ("bubbles"), leads to a 
linear system that is identical to the one obtained by using shape functions "^f , 
/ = 0, . . . ,A^, that satisfy ^f{xj) = djj and additionally 



( xj/fyv' - k^ xj/jv = Vv e //q' (k) n 

Since on a fixed mesh ,5^, we have limp^«, y/f = i//,, better stabiUty properties 
of higher order methods may reasonably be hoped for 
2. The system matrix of the Galerkin FEM based on the space V^'" is a tridiagonal 
matrix. The same matrix can also be obtained in different ways. Consider, for 
example, the sesquilinear form 

B{u,v):= / u'v' - kuv - iku{l)v{l) + Tk Lj^uL^v, (72) 

J Q. ^ J £2 
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where Lk = —-j^ — k^- For a suitable choice of the parameters Tk in dependence 
on k and Hk, the system matrix resulting from this B using the classical piece- 
wise linear hat functions leads to the same matrix as the Galerkin method based 
on the shape functions \//;, ; = 1 , . . . , A^. In ID, it is therefore possible to design 
nodally exact methods based on the stabilization techniques in the from (|72] i. 
In ifTSl . a nodally exact method is derived using other techniques. 

Acknowledgement: Financial support by the Vienna Science and Technology 
Fund (WWTF) is gratefully acknowledged. 



Appendix 

For the reference triangle K :~ {{x,y): < x < 1,0<3'< 1— -y} and /3 G [0,1) 
the following two lemmas require the spaces Hp^{K), Hp'^{K) as well as the 
Besov spaces B2^{K). The spaces B2^{K) are defined by interpolation using the 
/T-functional (see, e.g., ifTSl Chap. 12]). For m E {1,2}, the spaces Hp''"{K) are de- 
termined by the norm I1m||^;".;"(£) + ^^^^^ denotes 
the distance from the origin. 

Lemma A.l. Lef K be the reference triangle. Let j3 € [0,1). Then the embed- 
dings Hp'^{K) C B2 ^ (K) and H^'^ (K) C B2 ^ (K) are continuous. The embeddings 
Hp'^{K) CH--P-^{K) andHp'\K)cH^'P-^{K) are compact for all e > 0. 

Proof. Since the case j3 = corresponds to classical Sobolev spaces, we restrict 
our attention here to the situation j3 E (0, 1). The argument follows ideas presented 
in ifTTl Thm. 2.1] and ifTOI . We start with the following two Hardy inequalities for 
sufficiently smooth functions 

||r'^-^V«||^,(^) < C||«||^,.,(^), (A.1) 
||r'^-2(M-M(0))||^2(f) <C||«||^2.2(^); (A.2) 

here, jA.ll l is shown, for example, in |57, Lemma A. 1.7] and (IA.2b follows from 
combining ifTol Lemma 4.2] with jA.lb . Noting that ifTol (2.2)] states the continuous 
embedding Hp'^{K) C C{K), we have that «(0) in iA.2\ is indeed well-defined. 

We employ the real method of interpolation and write 00^ = {L^iH^)i-p /2.oo- 
Our method of proof consists in showing that for = 1— j8/2we have 

sup t^^K{t,u) <C\\u\\ 2.2,^-,, u:=u-u{0), 
te{QA) P 

for some C > independent of u. To that end, we proceed as in the proof of ifTOl 
Lemma 2.1]. For every 5 > 0, let Xs '= Cq{M^) with ;f = 1 on £5/2(0) and supp C 
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XdBsiO) as well as \\'^^Xs\\l'-{m.'^) < C5'', j e {0, 1,2}. We define the splitting 

Then from jAlT l and ( |A!2] i 

^ll«llL2((fnBg(0))\Bg/2(0))+<^^ 'll^"llL2((^nB5(0))\B5/,(0))+C|l^^«llL2(£\Ba/2(0)) 



From this, we can infer for any 5 E (0, 1) 

K{t,iT) < ||Ml||i2(^) +f|l«2||^2(£) < C\\u\\^2:. 



5^-P+t5-l' 



Selecting 5 = f'/^ gives K{t,u) < Ct^ ^^^\V^\\h--'^iPv Finally, the compactness 
assertions of the embeddings follows from the compactness of the embeddings 

^ioc C = for s' < s. □ 

Lemma A.2. Let j3 € [0, 1 /2) and K be the reference triangle. Then there exists C > 
such that for all u£Hp\K) there holds \\u\\^2f^gjl^ I^C llw|i/^2(^) + |j''^^M||^2(^) • 

Proof For each s> 1/2, we have the inequahty ||M||£2(5f) < CIImH^j^^^. From 
the embedding H^^^{K) C H^^^{K) of Lemma IaTI we then get ||«||i2(5^) < 

C\\u\\h~^(k) ^ C [||M|li2(jf) + \\r^^4L\K)\ ■ ° 

Lemma A.3. Let j3 G [0,1) and £2 dM^ be a finite sector with apex at the origin. 
Let M e C°°(i2) satisfy 

ll^n,/5,iV"+2M||^2(i2) < C„7„«! V« e No. 

Then, for k > ko > 0, there exist constants C, 7 > (depending only on j3, £2, Ju, 
and ko) such that 

ll^„,/5,iV''+2„||^2(i2) < CQ,k-<^^-P^fmax{n,k}"+^ V« G Nq. 

Proof. Lemma lA~4l vields 

max{«,fc}"+'' «! ^ 

where C, 7 > are independent of n and k. The result now follows. □ 
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Lemma A.4. Let /3 e [0, 1). Then for <r <R and all n G No 

\^ 'min{l,f±i}y max{n,k}"+^ - n"+^ 

Proof. We denote the left-hand side by Ihs and consider several cases. 
case 1: n <k and r{k + 1 ) < n + 1 : 



Ihs 



{k+l)r\"+^ 1 



« + 1 / 



for suitable C, 7 > if we assume that ^ > A-q > 0. 
case 2: n < k and r{k + 1 ) > n + 1 : 



111 1 



/t"+2 k^-P k^+P k^-P \ k J {k+lf+P 

< 



1 r 



Ihs 



k^-P \ k J \n+l 
<Cfk-^^-P^r"+P—^ 

for suitable C, 7 > 0. 
case 3: n > k: Then, for Q < r < R 



Ihs = (min{l,r})"+/^ -L- < r"+P -L- < k-^^-P)r"+P -L-n^-P 

^ L : J y „n+2 — „n+2 — ,,"+2 



<Ck-(^-P)r"+P^f 



for suitable C, 7 > 0. □ 



Further results and proofs 

Example A.5. In Example ll.il we studied the convergence behavior of the /i-FEM 
in the //I (12) -seminorm. In Fig.|4]we present the corresponding results for the con- 
vergence m the L^{a) -norm by plotting ||m ^ Miv||/^2(^)/||t'|lL2(^2^ vs. the number of 
degrees of freedom per wavelength A^;^ . For 79 = 1 , we observe 
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Il2 



which is in agreement with the analysis given in 1471 Sec. 4.6.4]. The cases p > I 
seem to behave differently as we observe 



\U-Un\\i2 



-(p+i) 



p=1 ; -u - \Cu = ^ 



p=2; -u - k' u = 1 
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Proof of Remark \4.9\ By interpolation using the /T-functional we can write 
^i+e ^ {H\H^)e.2 for e (0, 1). Hence, every u e H^+^{Q) can be decomposed 
as M = Ml + MT with 

ll«lll//l(I2) < f®ll«ll//l+e(i2)' Il«2||//2(i2) < t^^^hWH' + B^ny (A.3) 

where f > is arbitrary. The proof of Lemma |4~81 shows that ui and U2 can be 
approximated from Vn as follows: 
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^inf ll«2 - l'llL2(i2) V(m2 - V)||i2(i2) < /l^ll«2l|//2(i2) + {khf\\u2\\ i^2^i2) 

^inf ll«l -^'llL2(i2)+/^l|V(«l-^')llL2(i2) < /j|l«l|l//l(i2) + (^/0^ll«l|lL2(i2)- 



Using t =hm (IA.3l l we therefore get 

vg'vL " ^llL2(i2) V(m - V)||^2(i2) < /l' + ''||M||//l+e(i2) + (^/^)^ ll"lllL2(i2) + ll«2|lz.2(I2) 

The decomposition u = u\+U2 and the triangle inequality yield ||wi ||l2(^) + ||h2|Il2(^) < 
ll"llL2(i2)+2||Mi|lL2(i2) < ll«llL2(i2)+2|lMi|l^i(^) < ||M|1^2(i2)+2/i''|lM||^i+e(^). Com- 
bining these estimates, we obtain 

^nf Il"-V|li2(^)+/!||V(M2-V)||^2(i2) < {ll' + ^ + {khfh') ll«ll//I+e(^) + (Wi)2||M||^2(i2), 

which concludes the proof. □ 

Lemma A.6. Let j3 e [0, 1). Then, for every S N there exists a linear operator 
Tip : Hg'^{K) — >■ that admits an "element-by-element construction" in the sense 
of 1161] Def. 5.3] with the following approximation property: 

p\\u- 7tpu\\^2iK) + ll« - ^ Cp-^^-^'>\\rl^V^u\\^2^^y 

where C > is independent of p and u. 

Proof. Inspection of the proof of [61, Thm. B.4] shows that the operator Kp con- 
structed there does in fact not depend on the regularity parameter 5 > 1 . It has (as 
stated in |f6Tl Thm. B.4]), the approximation property 

pIIm- 7rpM||^2(^) + ||«- ^ C'/7"'''"''|M|ff-(£) e H\K), (A.4) 

if p > i — 1 . Upon writing the Besov space ^ as an interpolation space B2^ = 

{H^ {K),H^ (^) 1 for i G ( 1 , 2), we can infer for 5 e ( 1 , 2) from ( TOI l the slightly 
stronger statement 

p\\u- npu\\^2^K)'^h^ ^p"\\h^(k)^Cp'''''^^M\b\jk) Vm e (^) ■ ( A.5) 
Appeahng to Lemma lATI then yields 

p\\u-npu\\i2(^^ + ll« - < Cp'^'^^^Mh^^^ky ^^'^^ 

We replace the full Hj^ "{K) norm by the seminorm in the standard way by a com- 
pactness argument. Since Hp'^{K) is compactly embedded in //' (K) (see, e.g., ifTTl 
Lemma 4.19]) one obtains inf^^i^^ \\u — vW^i.i^^-j — C\\r^'V^u\\j^2(^gy The proof is 
completed by noting that (|A.4t implies that Ttp reproduces linear polynomials. □ 
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